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This  dissertation  describes  a procedure  for  arriving  at  a close  appnnsimation 
to  the  capadhuice  between  symmetrically  placed  conducting  strips,  possibly  of 
different  widihs,  on  opposite  sides  of  a dielectric  sheet.  The  procedure  is  based  on 
static  methods,  following  previous  work  by  Black  and  Higgins  for  total  capacitance 
of  the  structure  with  vacuum  everywhere,  and  Wheeler's  method  for  determining 


Dielectric  polaiiaation  is  included.  Refraction  at  the  vacuunt/dlelectric 


The  derived  equations  are  valid  for  all  finite  impedance,  all  values  of 


impedance  of  the  generalized  microatrip  on  a dielectric  sheet. 

The  methods  described  have  general  application  to  open  transmission  lines 
on  a dieiectric  sheet,  for  which  the  appropriate  conformal  transformations  can  be 
found. 

and  subroutines  for  the  basic  calculations  of  generalized  microstrip.  Tabulated  and 
graphical  results  are  induded.  Source  code  listings  of  the  demonstration  prograftis 


CHAPTER  I 
INTRODUCTION 


TEM  fTransverse  Electric  Magnetic)  transmission  lines.  TEH  transmission  lines 
have  a number  of  desirable  properties;  they  are  compact,  broadband  artd 
structurally  simple.  The  ground  plane  artd  dielectric  sheet  of  unbalanced  microsirip 
ate  assumed  to  be  of  infinite  width.  In  practical  terms,  they  are  simply  very  much 

analytical  microsirip  analysis  and.  in  particular,  analysis  of  microstrip  on  a 
dielectric  ^teet.  The  principal  transmission  line  parameter  is  characteristic 
impedance  (Z),  tvhile  the  structural  parameters  are  conductor  and  dielectric 
dimensions.  For  TEM  transmission  lines,  the  characleristic  impedance  is  fully 
determined  by  the  c^iacitance  per  unit  length  which  in  turn  is  obtained  from  the 
structural  paramclers.  At  this  writing,  the  analytical  difficulties  presented  by 
microstrip  are  still  far  from  resolved,  making  mictosirip  a priority  research  topic  in 
both  microwave  and  semiconductor  technologies.  It  is  the  accurate  delerminalion 
of  capacitance  from  microstrip  dimensions  that  presents  Ihe  major  obstacle.  The 
derivation  of  a close  analytic  approximation  for  the  capacitance  of  generalized 
microsirip  on  a dielectric  dteel  will  occupy  a major  pan  of  this  dissertation  and  is 
one  of  its  principal  contributions. 


Fig.  1-2.  The  asymmetrical  case,  shown  at  the  lop,  is  caJIcd  “unbalanced" 


microslrip.  The  symmetrical  case,  shown  at  the  boltom,  with  equal  conduciCFr 
widths,  is  called  'balanced'  microsuip.  Only  two  dimensions,  conductor  width 
and  conductor  spacing,  are  required  to  describe  these  cases.  In  fact,  either  case 
can  be  obtained  from  the  other  by  the  method  of  images.  Between  these  limiting 
cases  are  an  infinity  of  microstrip  configurations  where  both  conductors  are  of  finite 
but  unequal  width.  This  can  be  called  transitional  microstrip  since  it  provides  a 
tiansjeion  between  the  limiting  cases.  Description  of  transitional  microstrip  requires 
three  dimensions,  the  two  distinct  conductor  widths  and  their  separation. 

Discussion  of  generalized  microslrip  in  terms  of  three  dl.sciele  cases  is  for  purposes 
of  illustration  only:  the  Ihree  cases  are  totally  unified  by  the  derived  equolions 

Unbnlanced  microslrip  Is  the  most  common  and  universally  used  type  of 
miciosirip.  Many  devices  such  as  antennas  and  mixers,  however,  requite  a 
balanced  input  which  is  not  easily  oblained  in  unbalanced  miciosirip.  Generalized 
microstrip  provides  a connection  between  unbalanced  and  balanced  microslrip, 
enabling  direct  Uansforniation  between  the  two.  Such  transformation  structures  are 
called  balun  transformers  or  simply  haluns.  These  devices  occupy  an  important 
place  among  microwave  components.  Dcveli^raenl  of  precise  design  methods  for 

Microslrie  Analysis 

An  important  objective  of  microstrip  analysis  is  to  provide  accutale  design 
equations  to  support  a variety  of  microwave  and  semiconductor  dev 


vices.  Ideally 


this  might  consist  of  a ihree-dimensional  analytic  wave  solation,  yielding  the 
microstrip  dimensions  and  chaiacteristic  impedance  for  generalized  conductors  on  a 
dielectric  sheet.  Unfominately  this  ideal  corresponds  poorly  with  the  realilies  of 
microstrip  analysis.  Even  if  a formally  eaact,  analytic  wave  solution  could  be 
found,  and  it  probably  could  not  be,  it  would  almost  certainly  require  purely 
numerical  evaluation-  Numerical  wave  solutions  have  been  obtained  using  finite 
diflierence,  integral  equation  and  Fourier  methods.  These  methods  are  usually 
computationally  intensive  and  ultimately  require  comparisoa  to  an  analytic  soludon 
to  establish  their  validity  and  accuracy. 

It  is  an  olqcctive  of  this  dissertation  to  satisfy  as  many  of  the  chaiaclerisiics 
of  Ihe  ideal  solution  for  microstrip  as  possible.  An  analytic,  closed-form  solution  is 
sought,  since  its  accuracy  is  can  be  direedy  established  by  reference  to 
mathematical  and  physical  theorems. 


Microstrip  in  a lossless,  homogeneous  dielectric  with  lossless  conductors 
propagates  a purely  TEM  wave.  In  this  mode  there  is  no  axial  electric  or  magnetic 
field,  and  Ihe  two-dimensional  Laplace  (slatic)  solution,  In  Ihe  transverse  plane, 
satisfies  ihe  one-dimensional  wave  equation  fTelegrapher's  equation)  in  Ihe  axial 
direction.  For  microstrip  on  a dielectric  sheet  Ihe  didecuic  is  inhomogeneous,  with 
the  flux  partially  in  dielectric  and  partially  in  free-space.  In  general,  Ihe  flux 
impinging  on  this  boundary  is  distorted,  the  exception  being  flux  that  is  normally 


i-Chrisioffel 


the  dieleclric  boundary  conditions.  I^rictions  of  the  Schwan 
inuisformaiion  ut  sucb  that  the  gcneiaJ  boundary  conditions  cannot  be  imptsed;  a 
soluUon,  however,  can  be  obtained  by  neglecting  refiaction.  This  is  equivalent  to 
assuming  nortnaJ  incidence  for  all  flux  impingent  on  the  dielcctric/lree-space 
boundary.  This  is  the  only  analytic  assumption  required  to  obtain  a static  solution. 
Tlie  effect  of  this  assumption  is  assumed  to  be  small,  and  will  be  shovvn  to  be  so  by 
comparison  to  a numerical  solution  that  includes  the  eflecls  of  refiaction. 

The  solution  derived  in  this  disserution  consists  of  the  exact  solution  for  the 
static  capaciiance  with  vacuum  everywhere  and  a second  exact  soluUm  for  the 
poUriaaben  capaciiance  of  the  dielectric,  both  obtained  with  the  Schwarz- 
Christoffel  iransformaiicn.  The  mixed-dielectric  equations,  obiained  by  combining 
ihese  two  solutions,  are  analyiically  exact  except  for  refraction  cffecis. 

AxsuTnnrions 

Use  of  the  static  solubcn  as  a wave  solution  requires  ihe  assumption  thai  the 
microstrip  propagates  a purely  TEM  wave.  For  microslrip  in  a homogeneous 

dielectric,  however,  it  is  an  approximation.  This  well-known  approximation  is 
commonly  called  the  quasi-TEM  or  quasi-static  assumption  and  assumes  that  the 
axial  electric  and  magnebc  fields  are  negligible. 

Losses  in  the  conductors  and  dielectric  give  rise  to  an  axial  potential 
gradient  and  thus  an  axial  electric  field,  so  Ihe  conductors  and  dielectric  must  be 

ver,  copper  and  even  gold 


umed  to  be  lossless.  This  is  a good  assumption  for  silv 


conducton.  Similarly,  for  microwave  materials,  the  lossless  assumption  for  the 
dielectric  is  a good  one.  The  dietecuic  is  also  assumed  to  be  nonmagnetic,  a 

Use  of  the  Schwara-Christofrel  (conformal)  transformaiion  requires  two 
additional  assumptions:  (1)  the  effect  of  refraction  at  Ihe  dielectric/free-qnce 
boundary  is  assumed  lo  be  negligible  and  (2)  Ihe  conductors  are  assumed  to  be  so 
thin  as  to  approxlntate  malhemaiical  planes.  The  first  Is  necessary  because 
conformal  mapping  cannot  treat  general  dielectric  boundary  conditions.  The  sccortd 
is  to  simplify  Ihe  Schwarz-Chrisioffel  liansformalion  sufTicienlly  to  yield  a 
closed'fbrm  analytic  solution. 

The  quasi-TEM  solution  for  the  inhomogeneous  dielectric  case  does  not 
display  Ihe  frequency  dependence  (dispersion)  which  is  known  to  exist  In 
microstrip,  meaning  that  the  assumptions  are  imperfect.  They  are,  however,  quite 
good,  and  using  reasonable  design  nstrictions  so  as  not  10  sIgnificanUy  violate  the 
assumptions.  re.suits  have  been  excellent.  The  solution  is  exact  for  a,  = I and 
«,-><»  and,  with  Ihe  exception  of  refraction  error,  is  also  exact,  at  the  static  limit, 
for  all  other  values  of  relative  dielectric  constant. 

Related  Prior  An 

or  art  is  the  work  of  this  author  [t],  concerned 
electric  sheet.  The  homogeneous  dielectric  case 
ously  analyzed  by  Black  and  Higgins  [2],  using 


nformal  mapping. 


8 

A modiricadon  of  iheir  analysis  is  used  in  Lhis 
dissenalion  u>  determine  the  tola!  capacitance  of  generalized  microstrip  in  free- 

*nte  special  case  of  a conducting  strip  separated  from  an  infinite  conducting 
plane  by  a dielectric  sheet  was  analyzed  and  chatacterized  by  Dukes  [^-d]  with  the 
aid  of  an  tiectrolytic  tank,  and  later  by  Wheeler  [S-7]  using  approximate  conformal 
mapping  and  an  interpoiation  technique  which  produces  a correction  term  to  die 
resuit  obtained  by  Dukes  [3]. 

In  a 1984  paper.  Callarotli  and  Gallo  18]  describe  a method  for  analyris  of 

in  fiux-poiential  coordinates,  followed  by  a numerical  evaluation  of  microstrip 
capacitance  in  the  transform  plane. 

Dielectric  Green's  Function 

Weiss  er  at.  (9-10)  developed  computer  programs  treaUng.  as  one  cation, 
the  dielectric  Gmen's  funedon  as  described  by  Syivesler  [1 1],  Like  die  references 
included.  It  is  used  In  Chapter  4 to  esdmate  the  error  introduced  by  the  omission 


anaJydc 


Baiun  Transformers 


The  tapered  micrastrip  balun  ccndguiation  waa  described  by  Rumsey  (I^l- 
BaJuns  of  thia  type  can  simultaneously  provide  both  unbalanced  to  balanced  mode 
conversiwi  and  a suitable  impedance  transfonnalion  between  an  unbalanced  coaiial 
transmission  Ime  and  the  balanced  terminals  of  an  antenna  or  other  device.  The 
Klopfeostein  [I3|  impedance  taper  provides  an  c^timum  impedance  transformation 
but  requires  a detailed  knowledge  and  control  of  impedance  at  every  point  along  the 
taper.  The  procedures  provided  In  this  dissertation  ate  intended  to  accurately 
provide  those  parameters. 

This  dissertation  is  concerned  with  generalized  microstrip  of  the  type 
analyzed  by  Black  and  Higgins  [21,  but  on  a dielectric  sheet  as  treated,  for  qiecial 
cases,  by  Dukes  [3-4],  and  Wheeter  [5-7],  Enact  conformal  mapping  is  used,  as 
opposed  to  approsimate  conformal  mapping  empioyed  by  Wheeler  [5|. 
Disadvanu^es  of  this  approach  are  Ihe  added  complexity  of  the  exact  solution  and 

o^sei  by  Ihc  lack  of  appiosimation  in  the  resulting  solution. 

For  a given  characteristic  impedance  Z.  and  relalive  dielectric  constant  c, . a 
parameter  R determines  the  ratios  of  conductor  width  to  separation,  o/ft(R,Z.c,)  and 
bflKR.Z.eJ-  * = 0.0.  the  conductors  are  of  finite  and  equal  width,  wNIe  R 
« 1.0  gives  0 single,  finite  conducting  strip  over  an  infinite  ground  plane.  Values 
of  R between  these  limits  result  ta  two  finite  conducting  strips  of  unequal  widih. 
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II  would  be  desirable  lo  have  explicil  expressions  of  the  form 
R - R(aJh^.e,)  and  /!  = R(b/h,Z,e,l.  Il  ii,  however,  an  unfortunate  consequence 
of  the  solution  lhai  R arises  as  Uie  paiameier  of  an  elliptic  iniegral  of  the  third  kind 
and  cannot  be  expressed  in  this  way.  The  effect  of  this  implicit,  rather  than 
explicit,  reladceiship  between  conductor  width  and  the  parameter  /f  is  a somewhat 
more  difficult  calculation  than  would  be  ihe  case  if  an  explicii  relationship  could  be 


body  Is  devoted  primarily  to  development  of  the  physical  aspects  of  generalized 
microstrip  analysis  and  balun  design.  Mathemalical  developments,  such  as 
conformal  mapping,  are  treated  in  delail  in  the  Appendixes.  Source  listings  for 
computer  piograms  are  also  lisled  in  the  Appendixes.  This  organization  is  not 
intended  lo  diminish  Uic  imponance  of  either  the  conformal  derivalions  or  Ihe 
computer  programs,  but  rather  to  avoid  mixing  Ihc  physical  and  malhemabcal 

Of  general  inleresi  to  the  microwave  and  semiconductor  communities  are  the 
design  equntions  for  unbalanced  and  balanced  microslrip  in  their  most  accurate 
form.  With  an  estimated  accuracy  of  0.001  of  impedance,  these  equations  are 

Sons  provide  a standard  for  testing  numerical  microstrip 


thought  to  yield  ihe  me 


analysis  programs  and  assessing  their  accuracy  at  static  limits  as  well  as  for  direct 
use  In  the  design  of  unbalanced  and  balanced  microstrip. 


CHAPTER  2 

DESCRIPTION  OF  THE  CONFORMAL  TRANSFORMATIONS 
Backeroiinil 

Conromal  mi4)ping.  uid  the  Schwarz-ChrisiorTcl  transrormaiion  in 

such  as  Ihal  of  generalized  microsuHp,  inlo  a simpler  geamury  where  solutions  to 
the  Two-dimensional  Laplace  equation  can  be  more  easily  oblained,  usually  by 
inspeebon.  The  lesulbng  solutions  apply  lo  a uansverse  cross-section  of  (he 
structure,  with  the  microsirip  assumed  lo  be  consont  in  the  axial  direebon.  Direci 
applicabon  of  confonnal  mapping  is  restricted  to  cases  with  homogeneous 

microsirip:  bolh  are  for  homogeneous  dieleclric.  The  flrst,  derived  in  Apqxnrdix  A. 
is  used  to  obtain  the  total  capacitance  of  the  conduciing  strips  with  vacuum 

transformabon,  derived  in  Appendix  B,  is  used  lo  find  Uie  polaiizalion  capacimnee 
of  (he  dielectric  due  to  eleclrificabon  of  the  inner  faces  of  Ihe  conducting  sirips. 
I.e.,  the  faces  In  direci  contact  with  Ihe  dielectric  sheet.  This  is  also  a 
homogeneous  case  since  Ihe  polarizalion  field  is  eniircly  within  Ihc  dielectric. 
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The  Schwarz-ChrUiofTel  transrormation  canronDally  traiuronm  the  inlerior 

the  transformation.  In  space  coordinates.  Fi|.  A-l  of  Appendix  A,  the  Schwarz- 
Chrisloffel  polygon  is  Uie  upper  half  of  ihc  mlcrosirip;  it  encloses  Uie  conductors 

(j-plane)  polygon  thus  includes  infinity  as  one  of  its  vertices.  The  free-space  w to 
c-plane  transformation  is  obtained  using  two  setiueniial  conformnl  tmnsformaiions; 
the  first  transforms  the  r*plane  polygon  Into  Ihe  upper  half  of  [he  r*plane  while  the 
second  transforms  the  upper  half  of  Ihe  i-plane  into  Ihe  tv-plane.  or  flux-potential 

only  the  real  axis  of  the  r-plane  is  explicitly  used.  The  conductors  in  space 

outer  surfaces  face  the  inlerior.  The  Iransfonned  flux  lines,  in  the  iv-plane,  are 
lines  of  constani  v and  ihe  equipoiential  lines  are  those  of  constant  u.  When  the 
z-plane  is  mapped  inlo  a square,  as  in  Fig.  A-1,  the  squares  in  flux-poiential 
coordinates  map  bto  conformal  squares  in  space  coordinates.  This  is  shown  io 
Figs.  2-1  and  2-2.  where  each  square  in  Fig.  2-2  maps  to  a conformal  square  in 
Fig.  2- 1 . Presentation  of  Ihe  lu-plane  as  a square  is  done  to  illustrate  conformality 


nerally  repn 


+ 1.0 


i'-Axia,  X^x/b 

— Conductora, Plui  Unea, Equipotentlal  lines, 

Hux  division  line, Dielectric  boimdaries. 


Fi*.  2-1  GenersJixed  microstrlp  for  = 0.8.  end  Ireo-a 
impedance  108,4  Ohms,  space  coordinates. 


Cocductora, Flux  Unfls, Equipotentlal  lines, 

‘--Flux  diTleion  line,  — —Dielectric  bounderies. 


Fig.  2-2  Generelixed  micreatrip  for  A = 0.6,  end  free-epece 
impedence  166.4  Okuna,  flux-potentlel  coordinetee. 


Space  and  Flin-Poiaitkl  Coordinales 


While  it  ia  not  necessary  to  fully  understand  the  conformaJ  mapping 
derivations  of  Appeitdiies  A and  B to  follow  the  development  of  the  capacitance 
and  impedance  equations,  it  is  necessary  to  understand  the  nature  of  the 
transformed  (flux-potential)  parameters  and  their  relationship  to  generalized 
miCTOStrip  in  ^»ce  coordinates.  Figure  A-I  and  Table  A-1  in  Appendix  A show 
the  relationship  between  points  in  space  coordinates  and  their  images  in  flux- 
potential  coordinates-  These  relation^ips  ate  also  illustrated  in  Figs.  2-1  and  2-2, 
which  show  the  flux  lines,  equipotential  lines  and  dielectric  boundaries.  Parameter 
definitions  for  both  the  space  and  flux-potential  planes  and,  where  ^propriate,  the 
equation  number  for  calculation  of  the  parameter  are  given  in  the  Liat  of  Symbols. 

Figure  2-1  shows  Ihe  upper  half-plane  of  generalized  microslrip  with 
unequal  width  conductors.  The  space  coordinate  plane  is  normalized  by  Uie  half- 
separation  of  the  condtclors.  such  that  AT  = st/h  and  T = y/h.  The  conduciors  are 
the  bold  vertical  lines  alX  = -1.0  and  X = -H.O.  They  are  assumed  to  be  ol 
potentials  -1  V and  -H  V,  respectively.  The  thin  solid  lines  between  the 

upper  edge  of  the  narrow  conductor  and  on  the  inner  face  of  the  wide  strip,  is  the 


al  poCenlial.  The 
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e shaded  region  between  the  conductor]  indicates 
the  area  that  wiil  be  occupied  by  dielectric  but  is  currently  rilled  with  vacuum.  The 
venical  long-doubte-sbon  dashed  lines,  bounding  the  shaded  area  above  the 
conductors,  are  the  nanmeulliaed  surfaces  of  the  Inlinilely  wide  dielectric  sheet 
separating  the  conductors.  In  the  drawing  the  shaded  area  is  filled  (as  indicated)  by 

Flus-noiemial  coordinates 

Fig.  2-2  is  the  nus-polenUaJ  image  of  the  microsiiip  structure,  in  Fig.  2-1. 
resulting  from  the  conformal  transformation  derived  in  Appendia  A.  The  line 
coding  and  other  conventions  are  the  same  as  in  space  coordinates.  Three  points  in 
fliu-potential  coordinates  are  of  special  signilicance.  These  are  at  (u=^f,  v=o'). 
(B=-*r,  v=b’)  and  (u^Uj,  v=K').  Point  (/f,  o')  is  the  flux-potential  image  of 
point  1.0,  Y=a/h)  io  space  coordinates.  The  conductor  segment  at  u = *K, 
extending  from  v = 0-0  to  v = o',  is  the  flux-potential  image  of  the  inner  face  of 

point  (K.  o'),  and  extendmg  to  point  (AT.  1C),  is  the  flux-polcotial  image  of  the 
outer  face  of  the  narrow  conducting  strip  in  space  coordinates.  Point  (-K,  b')  has 
the  same  meaning  for  the  wide  conducting  strip.  Point  (Oo.  1C)  is  the  flux-potential 
image  of  Infinity  in  space  uoonlinatcs.  The  curved  long-double-sbort  dashed  lines, 

nonmelallised  surflices  of  the  infinite-width  dielectric  sheet  separating  the 


Legend: 

'■Conductore, ^Flur  Unee, Equipetentin]  linee, 

Dielectric  boundaries. 


Fig.  Geoeraiized  microstrip  for  R = 0.0  to  i.O.  and 

free-space  impedance  94.3  Ohms,  space  coordinates. 


Pic*  2~4  G^n^rftlizBd  ooicroBtrip  tot  It 
troo-sfoco  imped«noe  94.2  Ohms,  fli 
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flux-potential  coordinates  using  eqs.  (A-64)  and  (A-65)  for  the  general  case. 
Equations  (A-71)  and  (A-73)  can  be  used  for  the  special  case  of  balanced 
microsuip.  The  shaded  regions  in  boUi  figures  are  filled  with  vacuum  but  represent 
the  areas  that  will  be  filled  with  dielectric.  TTie  flux  and  equlpotential  lines  shown 
are  those  appropriate  to  free-space. 

The  solid  horizontal  lines  in  flux-potential  coordinates  are  flux  lines,  while 
■he  venicaJ  alternate  long-dton  dashed  lines  are  the  lines  of  equal  potential.  The 

inner  and  outer  faces.  Flux  below  this  line  is  from  the  inner  face,  while  that  above 
is  from  the  outer  face.  Limits  of  the  parameter  o'  are;  KH  & o'  & K.  Thus,  at 
least  half  of  the  electric  field  is  always  between  the  inner  faces  of  the  conductors. 
All  flux  from  the  outer  face  is  pariinlly  in  the  unshaded  (free-space)  region  and 
partially  in  the  shaded  (dielectric)  region  whereas  Dux  from  the  inner  face  is 


0.0  to  1.0  in  space  and  flux-potential  coordinates.  The  flux  and  equlpotential  lines 
have  been  omitted  in  Fig.  2-4  flux-potentinl  cooidinales  for  clarity  and  readability 
of  the  composite  graph.  All  six  illustrations  arc  at  the  same  impedance,  94.2  D. 
Both  conductors  stmulianeously  change  width  as  A Is  varied  to  irtaimain  this 

labeled  with  the  appropriate  values  of  I?  to  allow  correlation  with  their  space 


2-1.  Fi| 
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coordinsue  images.  The  lines  and  boandaries  are  coded  as  in  Fig.  2-1.  Figures  2-2 
and  2-4  Ulustraie  the  changes  in  the  Field  lines  and  dieleciric  boundaries  as  strip 
width  is  changed,  al  constant  impedance,  by  variation  of  R, 

Effect  of  the  Dielectric  Sheet 

The  introduction  of  dielectric  with  relative  permittivity  greater  than  one  into 
the  shaded  regions  has  two  effects.  The  more  important  effect  is  that  the  dielectric 
is  polarized  by  the  electric  field,  with  the  polarization  capacitance  rapidly  becoming 
the  dominant  capacitance  as  the  relative  dielectric  constant  is  increased.  The 
second  effect  is  refraction  at  the  vacuum/dielectric  boundary.  This  effect,  of 
necessity,  is  not  included  b the  calculations  but  Us  effect  is  estimated  b Chapter  4. 


AU  iUusliations  of  Reids  in  this  dissertation  are  either  vacuum  everywhere  or 
of  the  polarization  field  within  the  dielectric.  Bob  of  bese  are  homogeneous  cases. 

necessary  because  the  fields  cannot  be  ezpllcilly  calculaled  for  the  inhomogeneous 
case  by  be  mebods  of  bis  disseitalion. 

Figure  2-5  shows  the  polarizalion  field  and  equipoiential  lines  wibin  be 
dielectric  sheet  b be  upper  hnlf  plnnc  of  ^ce  coordinates.  The  relationships  bai 
delennine  bis  field  distribution  are  derived  in  Appendix  B,  The  case  illuslratcd  is 
for  balanced  microstrip.  The  figure  shows  twenty-one  field  lines,  equally  spaced: 
hence.  9S  percent  of  the  field  is  below  be  uppermost  field  line  and  only  5 percent 
lies  above  it.  The  dielectric  boundaries  are  also  the  oulermosl  field  lines;  bese 
extend  to  mfinity.  The  concentration  of  eleciric  field  al  the  conductor  edges 


ss 


Legend: 

Conductora, Flux  Uu 


Pig.  8-6  The  polarieetlon  Held 
conducUng  strips,  0.0,  is 


In  diele 
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for  relative  permittivity  gieater  Uian  one,  i.r„  only  when  a polarizable  dielectric  is 

For  a relative  permittivity  of  one  (free  space)  the  electric  field  terminatini 

in  Pig.  2-2,  and  with  no  polarization  field.  With  infinite  permittivity  the  electric 
field  terminating  on  the  inner  faces  of  both  conductors  is  that  shown  in  Fig.  2-S. 
Between  these  limits  the  field  on  the  inner  faces  is  a mishire  of  the  free-space  and 
polarizadfti  fields. 


Conformal  mapping  analysU  of  goieralizcd  microalrip,  and  in  particular  the 
Inhomogeneous  dielectric  case,  involves  specialized  lechniques  and  relationships 
lhai  are  not  often  encountered  elsewhere.  This  chapter  is  intended  to  introduce 
those  definitions,  relationships  and  techniques  required  fora  full  understanding  of 
the  derivations  unique  to  microslrip  analysis  by  conformal  mapping. 

MicrosiTin  Canaciunre  Calculations 

The  principal  reward  for  the  effort  involved  in  deriving  the  conformal 
transformalion  of  generalized  microstrip  is  the  ease  with  which  capacitrmce  can  be 

dielectric  constant  times  conductor  width  divided  by  conductor  separation.  This 
reladonship  is  well  known  and  widely  used  bul  in  ihe  inlerest  of  completeness  and, 
because  of  its  importance  within  this  dissertauon,  its  deiivalion  with  respect  to 
microstrip  is  included  here. 

The  w-plane  transformalion  of  microstrip  derived  in  Appendiz  A and 
gives  rise  to  an  electric  field  of  magnitude  £.  From  Gauss'  law  D = r£  = u and 


olLs/fneler 


(3-1) 


where  e is  Ihe  surface  charge  denslly  on  the  Inner  surface  of  one  of  the  conductors. 
Since  the  electric  field  U constant  in  magnitude  and  direction,  the  potential 
difference  between  the  conductors  is  given  by 

V,-V.  = £r,  Volts  (3-2) 

where  s is  the  sepanbon  of  the  conductors.  The  total  charge  per  meter,  Q.  on 

the  conductor  width.  Substitution  of  lE  for  u gives 

Q = iBw.  Coulombs/meler  (3-3) 

Capacitance  is  defirted  as  total  charge  on  either  conductor  divided  by  potential 
difference  of  the  conductors.  Dividing  (3-3)  by  (3-2)  gives 

e = Qnv,  - V,)  = e^,vts.  Farads/meier  (3-4) 

From  (3-4)  it  is  seen  that  capacitance  in  the  w-plane  is  calculated  as  dielecttiic 

«,  = I giving 

Co  = e^vils.  Farads  (3-5) 

Relative  Dielectric  Conslam, 

vacuum,  and  c when  a capacitor  with  the  same  geometric  capacitance  is  filled  with 
a dielecinc  substance,  the  relative  dielectric  constant  can  be  obtained  by  dividing 
(3-4)  by  (3-5)  giving 


(3-6) 


The  concqx  of  geometric  cspsciiance  is  based  on  an  Idea)  capacitor  whose 
free-space  capacitance  per  unit  length  is  completely  characieriaed  by  (3-S),  The 
divirion  of  (3-5)  by  (3-4),  which  lesults  in  (3-6),  can  be  viewed  as  a nomulizalion 


comaiiu  both  the  dimensioas  and  the  free-space  dielectric  constant,  the  result,  c„ 
can  be  interpreted  as  the  pseudo-capacitance  of  a unit  capacitor  with  w = s ^ ), 
whose  capacitance  is  dependent  only  on  This  interpretation  is  of  iitUe 
importartce  with  respect  to  relative  dieleciric  constanl  since  the  unit  pseudo- 
capacitor  is  urtifdmily  rilled  with  a single  dielectric  substance.  The  utility  of  the 
concept  is  in  the  calculation  of  effeclive  dieleciric  constant  where  the  unit  capacitor 
is  partially  rilled  with  a dieleciric  substance  and  partially  with  vacuum. 


If  an  ideal  capacitor  is  only  partially  rilled  with  a dielectric  substance  and 
partially  with  vacuum,  as  is  the  case  with  ihe  riui-polenilal  plane  representation  of 
microstrip  on  a dielectric  sheet,  the  definition  for  relative  dielectric  constant  can  still 
be  used  if  it  is  understood  to  represent  effective  dieleciric  constant.  Effective 
dielectric  conslanl  is.  therefore,  the  single  relative  dielectric  constant  rilling  Ihe 
entire  riua-potential  plane  that  produces  an  equivalent  capacitance.  Effective 


process  with  nonnalizaticm  constant 


:ince  this 


(3-7) 


effective  capaciunce  with  space  poniaJly  fillet]  wlL 


c„  is  the  Im-^ace  capacitance  of  the  microstrip.  For  application  of  Ibis  concept  to 
the  nut-poiemial  cepresenlatioii  of  microstrip,  the  flux-potentlaJ  plane  is  divided 
into  stibr^ons,  each  entirely  filled  with  either  dielectric  or  vacuimt,  and  with  their 
dielectricivacuum  boundaries  aligned  either  with  flus  or  equipolutuUl  lines.  This 

conformal  maps.  Each  of  these  subregions  is  then  normalized  to  and  each 
represents  an  apprc^riale  fractioD  of  a unit  capacitor,  with  the  sum  of  Ihe  areas  of 
the  subregions  equal  lo  unity.  Subregions  with  their  dielecuic/vacuum  boundaries 

equipotential  lines  are  in  series.  Each  subregion  is  a capacitor  with  homogeneous 
dielectric  whose  pseudo-capacitance  is  computed  as  conductor  width  over  conductor 
separation  times  either  e,  or  unity,  dqiending  on  whether  the  subregion  is  filled 
with  dielectric  or  vacuum.  The  effective  dielectric  constant  is  obtained  by 
approprialely  combining  the  pseudo-capacitances  of  the  subregions  in  series  or 

Eouivalent  Circuit  Comoonenls 

Derivatim  of  effective  dielectric  constant  follows  the  principles  set  down  in 
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curved  dieleelric  bmindehes  above  a'  do  nol  align  with  either  equipolential  lines 
(series  csmponem)  or  flui  lines  (parallel  cairiponenl)  but  are  a combination  of  both. 
Figuie  3*1  is  an  equivalent  circuit  of  microstrip,  showing  the  (our  capacitive 
components  that  must  be  determined  to  find  the  effeciive  dielectric  constant.  In 
Fig.  3-1  the  two  capacitors  in  series  represent  the  capacitance  attributable  to  flux 
originating  on  the  outer  face  of  the  narrower  strip.  This  will  be  called  the  external 
capacitance.  Conversely  capacitance  due  to  flux  originating  on  the  inner  face  will 
be  called  Ihe  inierrtal  capacitance.  The  concept  of  internal  and  exiemal 
capacitances  will  be  further  formalized  in  the  following  sections. 

Fractional  Parameters 

As  outlined  earlier  in  this  chapter,  the  effective  dielectric  constant  is  found 
using  fractional  pseudo-capacitances  of  subregions  obtained  by  normalization  of  the 
submgion  capacitance  by  the  total  free-space  capacitance.  For  consistency  it  will 
be  necessary  lo  define  fractional  width  parameters  as  well  as  the  fractional 
subregion  capacitance  parameters.  Fractional  widlb  parameters  are  defined  as 
partial  conductor  widlb  over  lout  conductor  width,  just  as  fractional  capacitance  is 
defined  as  partial  capacitance  divided  by  total  capacitance,  in  all  cases  the 
capacitances  will  be  'geonielric'  capacitance,  I.e.,  derived  as  conductor  width  over 
conductor  separation,  and  assumed  to  be  In  vacuum. 

Fractional  Width  Parameters 


2fl’/2S'  = 


(3-8) 


The  lefm  a’lK'  will  occur  rrequenlly  in  Ihe  derivalion,  so  lei 
A’  = a'llC. 

Similarly  the  fictional  width  of  ihe  ouler  face  is 


) - A'  = 2(K’  ■ fl')/2)r  = I - a'/K. 


(3-9) 


The  tola!  free-space  capacitance,  including  t>olh  upper  and  lower  half-planes, 
is  given  by  the  dielectric  ccnslani  of  free-space  times  conductor  width  over 


c,  = c^K'/2K  - i^lK. 

The  partial  parallel  free-space  capaciiance  of  the  inner  conducior  face  of  the  narrow 


6j2a'f2*:  = cjo'/s:.  (3-10) 

The  fractional  parallel  free-space  capaciiance  of  the  inner  lace  is 

A' = =a’IK\  (3-11) 

The  fractional  parallel  free-space  capacitance  of  Ihe  ouler  face  of  Ihe  narrow 


i - A'  = (e„2(/f  - o'niQlt^lJCi  = 1 - a'lK'.  (3-12) 

A consequence  of  normaliiation  is  that  the  fractional  widths,  capacilances 
and  areas  are  algebraically  equivalent.  They  are  nol.  however,  physically 
equivalent.  In  Ihe  derivation  which  follows,  the  physical  usage  of  Ihe  parameters 


Fractioial  Panllri  Didfaric  PnlthMHnn 


The  paiailel  dielectric  polariution  capacitance  is  derived  in  Appendix  B and 
given  by  (B-1 1)  in  nnnormalized  rorm,  derignaied  Cp.  For  use  as  a fractional 
capacitance  it  must  be  normalized  giving 

■ »-  (J-13) 

la 

CP  = QP<p)IP‘(p))HK‘IK).  (3-14) 


Cp/c,  = CF(e,  - 1).  (3-15) 

In  Fig.  3-1  this  is  the  parallel  dielectric  capacitance.  It  Is  analogous  to  the 
poraineler  designated  q'  by  Wheeler  [6].  Previous  to  Wheelers'  use.  Dukes  (3] 
used  the  parameter  implicitly  in  his  analysis  of  the  special  cate  of  equal  width 
strips.  The  value,  valid  for  the  special  case  of  equal  strips,  was  obtained  as  half 
the  capacitance  of  shielded  strip  transmission  line.  The  value  is  not  valid  for 
generalized  microsrrip,  but  is  exacUy  equivalent  with  R — 0.0  and,  by  images, 

R - 1 .0,  To  this  point,  two  of  the  four  capacitances  of  Ihe  equivalent  circuit 
have  been  found,  ^ifically  the  two  parallel  capacitances  compriang  Ihe  internal 
capacitance.  It  remains  to  find  the  series  capacilances  of  the  external  capacitance. 


Area  In  the  flux-poieniial  plane  has  little  meaning  when  taken  alone. 

Gearly  each  subregion  in  the  flux-potential  plane  has  an  area  and  Ihe  (olaJ  area  of 
these  suhr^ions  must  be  the  area  of  the  nux-potenlial  plane.  Beyond  this,  area 
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alone  has  Mule  meaning,  ll  is  Ihe  inlerprelalion  of  area  as  a capacitance  that  is 
meaningful.  The  capacilance  of  an  area  is  indeMnile  unless  either  width  or 
separation  is  spedMed.  In  the  derivation  of  effective  dielectric  constant  the  area  of 
the  dielectric  region  of  Ihe  flux-poleniial  plane  will  be  used  to  find  the  series 
dielectric  component. 

By  conservation  of  matter  Ihe  total  dielectric  area  in  Ihe  flux-potential  plane 
is  filed.  The  curved  dielectric  boundary,  see  Fig.  2-2,  gives  rise  to  both  a parallel 

directly  fbund.  The  parallel  component,  however,  is  known;  it  is  exactly  Cf*  given 
by  (3-14)  above.  1/  all  of  Ihe  dielectric  in  Fig.  2-2  is  assumed  to  be  in  parallel, 

parallel  dielectric  capacilance  which  will  be  designated  CPMAX.  It  is  the  maximum 
possible  value  of  parallel  dielectric  capociiance  with  all  dielectric  in  parallel. 
Meanwhile,  Ihe  actual  extent  to  which  the  dielectric  participates  in  the  parallel 
dielectric  capacitance  is  known  to  be  CP.  By  conservation  of  matter  Ihe  surplus 
parallel  dielectric,  CPMAX  - CP,  is  Ihe  dielectric  of  the  series  component. 

The  fractional  dielectric  capacilance  CPMAX  Is  analogous  to  the  parameter 

section,  the  designation  of  fractional  capacitance  is  thought  to  be  more  physically 
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Let  be  the  area  of  the  dielectric  in  the  upper  half-plane;  then  the  total 
dielectric  area  Is  2o,.  The  effective  width  of  Ihe  dieiectric  area  is  w = Tajs. 
where  s = conductor  s^iabon  - 2K,  giving  w - a^/K,  The  effective  paiaiiel 
free-space  capacitance  of  the  of  the  dielectric  region  is  e^ufs  — €qW/2K  — 

The  fractional  parallel  capacitance  of  the  total  dielectric  region  is 


The  area,  o^,  is  obtained  by  integrating  the  area  under  Ihe  curved  dielectric 
boundaries  in  Fig.  2-2.  This  gives 


where  v<u)  is  the  function  describing  the  dielectric.  The  integrand.  Ku).  is  given 
by  eq.  (A-64)  for  Ihe  general  case  and  (A-62)  for  the  case  It  = 0.0.  Use  of  (A-64) 
requires  on  iterative  procedure.  An  additional  difficulty  with  the  area  integration 
occurs  when  the  didecbic  boundary  becomes  double  valued  with  respect  to  a,  as 
shown  in  Fig,  2-4.  for  A 0.0.  When  this  occurs,  ^lecial  handling  is  required. 


The  total  microsirip  capacitance  can  be  represented  as  two  capacitances  in 
parallel,  designated  Q for  the  inietnal  capacitance  and  C,  for  the  external 
capacitance  (see  Fig.  3-1).  these  capacitances  are  defined  as  follows: 

C,  - the  capacitance  resulting  from  flux  terminating  on  the  inner  face  of  Ihe 


CPMM  = = aJ2KK’. 


(3-16) 


ninaling  on  the 


C,  - the  capaciunce  resulting  from  fluji  term 


e outer  face  of  the 


The  Imemal  capscilsuice  is  composed  of  two  fractional  capacitances  in 
parallel.  These  are  the  electric  component  and  the  dielectric  polarization 
component.  The  nature  of  these  capacitances  is  illustrated  by  reference  to  Bg.  2-2. 
The  internal  free-space  capacitance,  given  by  (3-11),  is  A'.  If  the  strips  are 


C.  = d'(e,  • 1)  + A'. 

The  first  term  on  the  right  is  the  dielectric  polarization  term  and  the  second  is  the 
electric,  or  free-space,  term.  For  the  homogeneous  case  the  polarization  and 
electric  fields  are  coirrcidenl  everywhere,  so  the  internal  capacitance  can, 
equivalently,  be  written  as  C|  w i/i'. 

The  polarization  field  of  Interest  in  microstrip  is  that  shown  in  Fig.  2-S.  For 
this  case,  the  polarization  field  is  constrained  by  the  finite  dielectric  boundaries. 

The  electric  and  polarization  fields  are  no  longer  coincident  and  the  Internal 
capacitance  must  be  eipressed  as  the  polarization  capacitance,  given  by  (3-lS),  and 
Ihe  free-space  capadlance  A',  giving 


C.  = CP((,  - 1) 


(3-17) 


Tliediff& 


CPMAX  - CP  represcnis  thai  portion  of  dielectric  that  does 


not  participate  in  the  irtie  parallel  capacitance  derived  in  Ai^endix  B and. 
therefore,  is  assumed  to  be  the  series  dielectric  component.  The  difference 
represents  an  estra  conductor  width  in  normalized  coordinates  with  a conductor 
separation  of  unity.  It  can  therefore  be  directly  interpreted  as  a dielectric  area. 

Figure  3-2  is  a schematic  representation  of  the  dielectric  and  free-t^ace 
areas  comprising  the  internal  and  external  capacitances  in  parallel.  The  external 

width  I -d’,  at  the  top  of  the  schematic.  The  vridlh  is  that  given  by  (3-9).  The 
traded  legion  on  the  upper  left,  of  separation  S and  width  W = \ • A',  represents 
the  series  dielectric  component  of  the  external  capacitance.  The  unshaded  region 
with  separation  I - 3 is  the  series  fiee-space  capacitance.  With  (he  width  of  (he 
series  dielectric  eqiial  to  I - ri'  and  its  area.  CPMAX  ■ CP,  the  separation.  S,  can  be 
obtained  from 

WXS  = CPMAX  ■ CP, 

giving 

S = (CPMAX  ■ CP)IW.  (3-lg) 

The  inverse  of  the  external  capacitance  is 


]/(i,WISl  + l/(H'f(l  -Xf) 


(3-19) 


Pig.  3-2  Normalized  w-plaoe  ehowiQg  fractional  capaoitancea. 


(3-20) 


yields  the  eipiesaon  for 
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' for  W and  the  right-hand  side  of  (3-18)  for  S in  (3-19) 


1 (,  \^CPMAX-CP  ' (3-21) 

maximum  possible  value  of  exlemal  capacicance. 

e.--  ] CPMAX-  CP  ' (3-22) 

Effeclive  Dieleclric  Conslanl 

uitemal  capacilance.  (3-17),  and  Ihe  exlemal  capacilance,  (3-21),  in  parallel  giving 
^ ^ CPMAX-  CP  ' (3-23) 

For  calculation,  subaituie  a'lK'  for  A‘,  where  s'  Is  given  by  (A-97)  for  the 
general  case  and  (A-104)  for  die  case  K = 0.0.  Compute  Causing  (3-14)  and 
CPMAX  wiOi  (3-16). 

It  is  noteworthy  that  the  width  of  Ihe  wide  conducting  strip  does  not  appear 
explicitly  in  the  equation  for  It  does,  however,  influence  Ihe  value  implicitly 
through  via  the  calculations  for  CPMAX,  CP. 

Three  derived  paramours  are  Involved  in  the  calculadon  of  1,,^  Tliey  are: 
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A‘(K,Z^  = ihe  rracdoiul  parallel  free- space  capacitance, 

CPMA)^S,Z^)  = Ihe  maximum  value  of  fractional  parallel  dielectric 
capacitance  and 

CnRAH  = the  true  value  of  fractional  paralld  dielectric  capacitance. 

Their  behavior  aa  timciioiis  of  free-apace  impedance,  Z^.  and  Ihe  parameter 
It  ia  shown  in  Figs.  3-3,  3-4  and  3-S,  respectively. 


The' 


is  completely 


determined  by  Ihe  eff 
that  the  dielectric  she 
equal  to  Ihe  permeabi 
microwave  dielectric 

The  characteristic  impedance  (Z) 


This  is  nearly  always  t 


'-ii- 


: characteristic  impedance  of  a microstrip  In 


where  p,  = arXlO"’  Henry/meler,  Ihe  permeability  of  free  space,  and  Cq  is  Ihe 
geometric  capacitance  per  meter.  In  Ihe  presence  of  a nnn-magnelic  dielectric  sheet 
between  the  conductors  the  characteristic  impedance  can  be  written 


(3-25) 


Pig.  3-3  NormBiiaed  inleraal  free-spacs  capacitance  (A‘)  Tcraus 
free— apace  Impedance  and  the  parameter  R, 


Pi*.  3-4  Nom4ll»ed  nmiimiim  parallel  dielectric  capacitance 
(CPMAX)  vereua  free-epace  impedance  and  the  parameter  R. 


FI*.  3-6  NornMiUz«d  true  pereUel  dielectric  capeeltoQoe  (CP) 
feee-apBce  impedance  end  the  peremeter  R. 


Dividing  (3-25)  by  (3-24)  cancels  giving 


(3-26) 

where  Cq  is  Ihe  geometric  free-space  capacitance  of  the  microstrip,  given  by  Eq. 
(A-li)  of  Appendix  A,  ^ is  the  free-space  characteristic  impedance  given  by  Eq. 
(A-I2),  c^  is  the  effective  microstrip  capacitance  on  a dielectric  sheet,  and  Z is  Ihe 
characteristic  impedance  of  the  microstrip  in  the  ptesence  of  the  dielectric  sheet. 
Since  = c^Co,  (3-26)  can  be  written 


The  most  complex  and  Bme  consuming  procedure  in  the  use  of  (3-23)  is  Ihe 
computation  of  CFMAX,  by  numerical  integration.  Wheeler  (6)  obtained  this  value 
by  approximating  Ihe  dielectric  boundary  vrilh  an  ellipse.  This  method  gives  good 
accuracy  in  the  mid-range  if  impedance  for  R = 0.0.  Because  of  the  wide  range  of 
boundary  shapes  occurring  in  generaliacd  microstrip  (see  Fig.  2-4),  this  procedure 
cannot  be  used  in  general.  The  following  emjarical  formula  docs  not  require 
integration  and  leads  to  Ihe  calculation  of  Ihe  characteristic  impedance  of 
generaliied  microstrip  on  a dielectric  sheet  within  less  than  ± 0.25  percent  of  the 
value  obtained  using  (3-23): 


<3-28) 
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1 - 0-536(1 


EqiiaUon  (3-27)  is  valid  lo  the  speciOed  accuracy  for  0 < Z,  □.  0 S 1?  5 I 
and  I S e,.  It  is  the  empirically  adjusted  average  of  the  cases  vrith  no  dielecinc 


CHAPTER  4 

COMPARISON  TO  PREVIOUS  RESULTS 


The  analyses  of  Dukes  [3]  and 
rroslrip  hul  for  the  ^eciaJ  eases  R 


are  not  for  generalized 
: = 1,0.  Their  fomulalion  of 

To  compare  the  equaUoiu  derived  by  Dukes  and  Wheeler,  the  eziemal 
capaatance,  eq.  (3-21),  can  be  written  in  the  foim 

(*-1) 


Since  X is  always  less  than  I.  (4-1)  can  be  expanded  in  a binomial  sc 


C,  = (1 -A')(i  + J(  + J? + ..,),  (4-2) 

Use  effective  dielectric  constant  can  be  approximated  using  one  or  more 
terras  of  (4-2)  as  the  external  capacitance,  in  parallel  with  the  intemaJ  capacitance, 
given  by  eq.  (3-17). 
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Cornnarison  wiUi  Dukes 

Uting  only  the  first  term  of  (4-2)  gives 

e,„  - CP{t,-l)  + 4’  -I-  1 - CP(«,-  1)  -I-  1.  (4-3) 

Eicept  fer  notation,  this  is  precisely  the  equation  derived  by  Dultes  [3], 
Effectively,  this  formulation  simply  ignores  the  series  dielectric  portion  of  the 
estemal  capacitance.  Z computed  using  (4-3),  and  using  the  exact  conformal 
equations  from  the  appendixes,  is  accurate  to  better  than  0.75  percent  of  Z, 
computed  with  (3-Z3).  This  shows  that,  at  least  for  the  balanced  mlcrostrip  case, 
the  dielectric  in  the  external  region  is  not  very  important. 

Including  the  first  two  terms  of  (4-2)  gives 

Crf,  - CP(4,-  1)  -h  4’  -h  1 -4'  + (1  - Ut,)(CPMAX-  CP). 

= CPU,  - 1)  -h  1 + (1  - ]I(,)(CPMM  - CP).  (4-4) 

Wheeler's  (6]  expression  for  consists  of  two  equations.  The  first  is  an 
interpolation  to  obtain  the  effective  'filling  fraction,'  and  the  second  to  obtain  the 
effective  dielectric  constant.  These  equations  are 


‘.e  = - 1)  + 1.  (4-6) 

Substitution  of  q from  (4-5)  into  (4-6)  gives 


g’U,  - 1)  + I + (I  - lU,)(q’  - q’). 


(4-7) 


In  Wlie 


ecler's  nouiion,  g'  U analogous  lo  CPMAX  and  q’  is  analogous  10 
CP.  WiOi  this  notsUonal  change,  (4-7)  is  eaacUy  equivalent  lo  (4-4),  showing 

solution  derived  in  this  paper  for  the  case  R = 0.0. 

Wheeler  dcnved  his  equations  using  approximate  conformal  mapping  rather 
than  exact  conformal  mapping  as  is  done  in  this  dissertation.  Numerical  tesults, 

Z,  compared  to  (3-23).  When  (4-7)  is  evaluated,  however,  with  the  exact 

hundredths  of  a percent  of  the  result  from  p-23>. 

Both  Dukes'  and  Wheeler's  solutions  are  for  the  special  case  R = 0.0. 
balanced  microstrip,  whereas  the  solution  derived  in  this  paper  is  for  geneialized 
microstrip.  The  solutions  are  conceptually  analogous  with  respect  to  the  series 
expansion  but  are  by  no  mans  equivalenL 

In  addition  lo  the  analytical  comparison  with  previous  static  solutions,  it  is 
instructive  to  compare  numerical  results.  A complete  comparison  for  the  general 
case  of  microstrip  is  not  possible  because  no  solution,  other  than  that  derived  here, 
exists.  The  case  R - 0.0.  balanced  microstrip,  and  its  image  countcipart, 
unbalanced  microstrip,  with  «=  1 .0,  can  be  compared  since  a number  of  solutions 
fOr  these  cases  exist.  Actually,  the  underlying  concepts  of  both  Dukes'  and 
Wheeler's  analyses  can  be  directly  applied  lo  the  general  case  using  the  general 


study  has 
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mapping  uquatuxis  derived  in  this  dissenation.  Such  a 
been  formally  done. 

Three  specific  methods  will  be  compared.  They  are;  I)  Dukes'  foimulaDon 
[3|  2)  Wheeler's  ^proximate  conformal  formulation  [7]  and  3)  the  dielectric 
Green's  function  formulation  of  Weiss  er  of.  [101. 

Dukes'  formulation  is  included  because  it  shows  the  relative  imporunce  of 
the  series  component  of  dielectric.  Recall  that  Dukes  omits  the  series  dielectric 
component.  Comparison  to  Wheeler's  approximate  conformal  method  illustrates  the 
dilTerence  between  approximate  and  exact  conformal  methods.  The  dielectric 
Green's  function  foimulntion  is  of  interest  for  several  reasons.  It  is  based  on  sound 
analytic  principles,  it  includes  refraction  effecht  and  finally,  the  computer  program 
is  available. 

In  the  comparison  with  Dukes  and  Wheeler,  which  are  both  conformal, 
differeiices  are  due  to  tyiproxiiTialions  made  by  the  analysis  while  difTetences  from 
the  dieleciric  Green's  function  result  from  numerical  and  discreiizatioD  errors  in 
computation  and,  to  a lesser  extent,  from  omission  of  refraction  in  the  exact 
conformal  mapinng  solution  of  this  paper.  The  dielectric  Green’s  function  method 
is  the  most  important  of  the  three  rince  it  includes  refraction.  A modified  form  of 

be  used  later  to  estimate  the  error  in  the  conformal  solution  resulting  from  omission 
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Figure  4-1  is  a direct  comparison  of  Z.  compulcd  using  <3-23)  of  this  paper 
with  Dukes’  formulation  as  given  by  cq.  <4-3).  Figure  4-2  compares  eq.  (9)  of 
Wheeler  [TJ.  and  Fig,  4-3  Uic  compuler  program  MSTRIP2,  of  Weiss  et  a/-.  [10], 
Both  Wheeler  [7]  and  MSTRIP2  differ,  at  most,  by  slighUy  more  lhan  1 percent. 
The  maaimuni  difference  (for  generalized  microstrip)  from  Dukes  (3)  is  slightly  less 
than  0.75  percent. 

Identification  of  the  difference  from  Weiss  er  at.  |I0]  as  numerical  or 
discretization  errors,  rather  than  errors  due  to  omission  of  refraction,  is  dreie  by 
comparison  at  c,  - 1.0  where  the  conformal  solution  is  enact. 

Estimate  of  Reftaeiion  Error 

The  totat  difference  between  the  presence  and  absence  of  dielectric  in  the 
eatetnal  region  is  small,  so  the  effect  on  impedance  of  refraction  at  the  dielectric 
boundaries  must  be  much  smaller  yet.  To  quantify  the  error  introduced  by 
omission  of  refraction,  a second  solution,  of  equivalent  precision  to  the  conformal 
solution  but  including  refraction,  is  needed.  The  difference  in  the  results  of  these 
solutions  would  then  appmnimate  Ihe  refraction  error  in  the  conformal  solutim. 

The  dielectric  Green's  function,  as  implemented  in  MSTRIF2,  includes 
refraction  but  numerical  and  discretiation  errors  are  so  large  as  to  obscure  its 

conformal  solution  of  this  paper  in  Ihe  range  e,  = 1 to  e,  -•  <»,  the  dielectric 
Green's  function  can  be  asymptotically  corrected.  Neither  the  conformal  solution 


Fig.  4-1  Uloroatrip  Impeduce  computed  b;  the  method  of 
Dukes,  oq.  (4-3),  end  by  eq.  (3-23)  of  this  dieeerUUoa 


Cofi/omul  lUp 


»/h 

Fi*.  4-8  Micro*trip  impeduce  computed  by  the  epproximete 
method  of  Wheeler,  eq.  (9)  of  Bof.  7,  end  eq.  (3-83)  of  thle 
dieierletioo  Tereua  conductor  width  end  rel,  diel.  conet. 
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CoDlonaiAl  Xap 


4-3  Mlcroitrip  Impeducs  cookpuCsd  by  tie  method  of 
U3TREP2  end  by  eq.  (3-23)  of  thje  diesertetion  ooreuo 
conductor  width  end  reUtive  dielectric  oonetent. 
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of  this  paper  nor  MSTRIP2  Impose  a limilalion  on  ihe  upper  value  of  dieleclric 
consont,  so  any  very  large  value  of  e,  within  the  noating  point  range  of  the 
computer  can  be  used. 

The  solution  of  this  dissertation  provides  the  enact  values  of  and  a//i  al 
e,  = 1.  The  solution  is  also  exact  for  very  large  At «,  ” lO"*  the  error  is  on 
the  order  of  Kr‘“  which  cenainly  can  be  considered  negligible. 

Procedure  for  Calibratine  MSTRIP2 

The  following  exact  parameters  from  the  conformal  solution  will  be  used: 

I.  & tf  ^ I,  Zq  (exact),  a/fi  (exact). 

II.  9 c,  “ 10'°.  Z - (exact),  a/h  (exact). 

1 . With  inputs  of  ^ and  with  e,  1 , find  alh  with  Ihe  conformal  procedure. 

NCTTE:  To  oblain  values  of  a/h  suitable  for  precise  comparison  with 
MSTRIP2,  the  value  of  free- space  impedance  used  in  Ihe  conformal  procedure, 
nominally  376.7  (1,  was  set  to  correspond  to  ihe  vdocity  of  light  used  in  MSTR1P2, 
2.99792458  X lO"  m/%. 

2.  With  inputsof  £,  = l,and  = a/A  to  MSTRIP2,  find  Z.  Let  this  be  2;^. 

3.  With  inpuu  of  2^,  and  e,  = lO'",  find  Z with  the  conformal  procedure. 

4.  Wilh  inputs  of  e,  = I0'“  and  w/A  = a/h  to  MSTRIP2,  find  Z.  Let  Ihis  be 
Z. 

5.  Form  Ihe  ratios  r’  = and  r"  = Z/Z..  (4-8) 

6.  Find  r for  interraediale  values  of  c,  by  interpolation  as  1/e,  using 


(4-9) 
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7.  With  the  value  of  alh  used  ahove,  and  any  value  of  I < c,  < 10'°  as  input 
to  MSTRIP2,  compute  the  value  of  Z,.  Adjust  Uie  value  of  Z,  using 


3.  With  the  value  of  Zg  used  alsove,  and  the  value  of  1 < c,  < I0'°  used  in  the 

ae,)- 

9.  Compute  the  percent  diffeience,  C,  as 


The  value  of  Z.(adjusied|  from  MSTRIP2  is  always  slightly  lower  than  the 
value  of  Z computed  by  the  conformal  method.  The  maximum  diffeience  of  d = 
0.0374  percent  occurs  at  or  near  a = 111.0  Q<o/A  = 1.3183874)  and  e,  = 3.6. 
This  difference  is  tentatively  attributed  to  the  omission  of  refraction  effects  in  the 
conformal  computation  ofZ. 

10.0  s Zg  S SOO.O.  (1  using  MSTRIP2  with  m = 20.  MSTRIP2  does  not  support 
conducting  strip  over  an  infinite  ground  plane,  K = 1.0,  are  thought  to  be 


Z.(adjusted)  = rZ,. 


(*-10) 


(4-11) 
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Edwards  [14]  gives  formulas  for  the  cutoff  frequencies  for  die  TM  and 
transverse  resonance  inodes.  Edwards  cites  these  formulas  from  Vendeiin 
(reference  37  of  Edwards  book)  who,  in  turn,  attributes  Collin.  The  formula  for 
TM  cutoff  Is 


where  c is  the  free-qace  velocity  of  light  and  r ::  2A, 


(4-12) 

thickness, 


Transverse  resonance  cutoff  is  opprosimatcd  by 

/cm* — = 

(2h..0-8i)^ 

with  w - 2o- 


(4-13) 


Relaiionshin  to  the  Wave  Eouation 

Assuming  that  the  condidoni  implied  by  (4-12)  and  (4-13)  are  met,  the 


simple  assumptions.  The  primary  source  c 
to  the  wave  solution  is  the  inhomogeneous 

of  the  conductors,  (f  the  series  component 
propagation  is  at  the  velocity  of  light  in  fre 
terminating  on  the  inner  faces  is  entirely  in 


1 the  static  solution  with  respect 

ng  on  the  inner  and  outer  faces 
sternal  field  is  neglected. 
Conversely,  the  ffeld 
c with  a velocity  of  propagation 
lists  independent  of  the  other. 


(4-14) 


: dielectric.  Using  (4-14)  the 


rcinlive  velocity  difTerEntial  between  the  inner  and  outer  fields  can  be  defined  as: 
i(e,)  = (c  - Vj)/c  - 1 - vp.  (4-15) 

The  reUtive  velocity  differential  is  seen  to  be  zero  for  e,  = 1 and  approaches  a 
maximum  of  1 as  e,  is  increased  to  a very  large  value.  The  iniemal  and  external 
capacitances  are  given  by  eqs.  (3-17)  and  (3-21),  respectively,  at)d  their  ratio  is  a 
measure  of  the  Held  fractions  terminating  on  the  inner  and  outer  faces  of  the 
conductois.  Define 

f(0  = CJC;.  (4-16) 

From  e().  (3-22)  the  external  capacitance  (CJ  has  a finite  maximum  value,  which  it 

limit  as  c,  increases  without  limit.  As  e, -•  o>,  f(e,)  -«0,  meaning  that  the 
external  field  is  negligible  with  respect  to  the  internal  field  for  very  large  {,. 

At  the  limit  for  low  impedance  the  conductor  widths  increase  without  limit, 
ultimately  becoming  an  infinite  radial  transmission  line  with  the  field  entirely  on  its 
inner  faces.  Therefore,  at  low  impedance,  the  velocity  differential  is  meaningless 
since  them  is  no  freld  on  the  outer  faces,  if  it  is  assumed  that  the  error  in  the  static 
approximation  is  proportional  to  8(«,)  and  to  the  ratio  of  the  external  U>  internal 
capacitance  f(«,),  the  estimate  of  errtrr  in  the  static  approximation  with  respect  to 

eiior  = c,S(«,)f(«,)  (4-17) 

wheie  C|  is  an  undetermined  proportionality  constant. 
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Calculations  made  with  (4-18)  give  Uie  error  relative  to  Hie  unknown 
conuani  r,,  To  Rnd  the  absolute  error  it  is  necessary  to  evaluate  the 
proportionality  constant  which  requires  a wave  solution  not  currently  available. 
Evaluation  of  (4-18)  for  the  relative  error  is  of  some  interest  since  it  indicates  the 
lelative  magnitude  of  error  as  a function  of  ZO.  /)  and  e,.  Also  since  c,  is  a 
particular  single  constant,  the  location  of  maiimum  of  the  relative  error  will  also  be 
that  of  maximum  of  the  absolute  error.  Numerical  results  for  these  calculations  will 


CHAPTER  5 

GENERAUZED  MICROSTRIP  CALCULATIONS  AND  RESULTS 

In  the  preceding  Chapters  and  the  Appendixes  A and  B.  about  two  hundred 
equations  have  been  derived.  The  thirty  six  most  Frequently  employed  are 
collected,  by  category,  in  Table  5-1  For  reference.  Two  demonstration  programs 
used  in  preparation  of  this  dissertation,  entitled  MISTI  and  M1ST2,  together  wilh 
their  supporting  subroutines  and  function  subroutines,  are  listed  in  Appendix  C. 
These  programs  are  written  in  double  precision  Fortran  77  and  are  designed  for  use 
wilh  “G"  floating  point  format.  They  should  also  function  in  “D"  floating  point 
formal  but  possibly  over  a slighUy  reduced  range  of  impedance.  Tbe  programs  and 
subprograms  are  extensively  commented  and  are  self  documenting. 

Proeram  MISTI 

Progmm  MISTI  .FOR  accepts  the  free-space  impedance  Zg,  the  parameter  X 
and  the  relative  dielectric  constant  e,  as  inputs  and  returns  the  coire^nding  values 
of  d/A.  b/A,  Z.  vp  and  the  transmission  line  capacitance  per  meter  in  PF  for 
generalized  microstrip  on  a dielectric  sheet.  The  calculation  flow  is  stiaightforward 
and  is  readily  followed  from  the  program. 
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TABLE  5-1 

SUMMARY  OF  FORMULAS 


TABLE  5-1  (Condmisd) 


wiihlAT'/I: 


I [ft’  ea-l. 


(A-58) 


v-IS- 


Formulm  for  Ihe  Sotgial  Cue  S - 0 


Traniform  Funelio 


(A-63) 


TABLE  5-1  (Continued) 


' Traiuforni . 


. F{T..k-). 


To  find  u ■ ulY.vJtl 


. I 


‘■‘litao-i'-flV). 


(A-71) 

(A-70) 

(A-72) 

(A-7J) 

(A-75) 

(A-74) 

(A-76) 


(A-102) 

(A-103) 

<A-I04) 


TABLE  S-I  (Continued) 


Difitciric  Polifitmlion  Ci 


c -e,^(e  -1) 


/ " ^sach^t  + accb’fi  - aechM  tect^B , 
where  ^ - ni/4A.  end  A - e«/4A. 

Effective  Dielectric  Cointatil 


(A-105) 


i-aswt-— 

Chetmclerriric  tmnedince  on  i Dielectric  Sheet 


(3-2«) 


(3-27) 
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Figure  S-l  is  a graph  of  generalized  microstrip  Idr  values  obtained  for  Z 
from  10.0  to  200.0  Ohms  and  it  from  0.0  to  0.99  with  e,  = 9.6,  99.S  per  cent 
alumina.  Any  number  of  such  graphs  for  general  values  of  A,  and  various  values 
of  dielectric  conslani,  can  be  computed  with  high  accuracy. 

Tables  5-2  and  5-3  are  estimales  of  the  maximum  relatiYC  error  of  Ihe  static 
solution  with  respect  to  a wave  solution.  The  values  of  maximum  relative  enor, 
given  in  Table  S-2,  are  relative  to  an  undetermined  constant.  Table  3-3  indicates 
lhai  Uic  maximum  difference  between  the  results  of  this  di.ssenation  and  a wave 
soluiion  should  occur  near  a relative  dielectric  constant  of  2.3.  This  information  Is 
useful  as  a beginning  point  for  future  comparison  to  a wave  solution  and  also  to 
avoid,  wherever  possible,  the  region  of  potential  error. 

Ptc^ram  MISTI  computes  Ihe  static  parameters  of  generaUzed  microslrip  in 

CPMM.  In  addition,  the  program  umullaneously  computes  impedance  using  the 
approximate  formulation,  eq.{3-28),  and  the  percentage  difference  between  the 
resulting  impedance  values.  The  peak  difference  is  shown  to  be  less  than  0.25  per 
cent,  more  than  adequate  for  all  practical  calculations.  The  computationally 

eq.(3-28).  This  can  result  in  a signiricanl  simpliricalian  and  speed  Increase  for 
generalized  microslrip  design  programs. 


63 


{sneriJiaed  microetrip  oa  &lumioa. 


TABLE  3-2 

MAXIMUM  RELATIVE  ERROR  VERSUS  ZO  AND  R 


20 


420.000 

210.000 

105-000 


26.230 

13.125 


0.145  0.145  0.145 
0.136  0.137  0.139 
0.103  0.110  0.U7 


0.145 

0.140 

0.125 


0.147 

0.145 

0.134 


0.036 

0.031 

0.017 
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TABLE  5-3 

(,  FOR  MAXIMUM  RELATIVE  ERROR  VERSUS  ZD  AND  R 


ZD 

840.000 

420.000 

210.000 
lOS.OOO 
52.500 
26.250 
13.125 


2.227 


2.227 

2.212 

2.227 

2.263 

2.292 

2-312 


2.227 

2.207 

2.212 

2.260 

2.309 

2.339 


2.226 

2.203 

2.207 

2.269 

2.321 

2.339 

2.334 


2.226  2.239 
2.203  2.227 
2.216  2.220 


2-275  2.251 
2-312  2.285 


2.317 


2-304 


0.6 
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Program  MIST2.F0R  accepts  the  desired  impedance  (Z)  of  the  microstrip  on 

conductor  widlii  alh,  vp  and  the  transmission  iine  capacitance  per  meter  in  pico 
Farads,  seleciably  for  either  unbalanced  or  balanced  microstrip.  Since  the  input  is 
the  impedance  on  a dielectric  sheet  2 rather  than  the  program  iterates  to  the 
appropriate  value  ofZa,i  method  that  is  required  for  many  derign  calculations  in 
generalised  microstrip.  As  was  the  case  for  MISTI,  this  program  implements  the 
most  accurate  form  of  the  derived  equations.  Integration  of  dielectric  area  is  by  use 
of  the  inverse  ttansformation  functions,  eqs.(A-TO),  (A-71),  (A-72)  and  (A-73),  for 
finding  the  dielectric  boundary.  These  equations  do  not  requite  iteration  and  are 
suffictenily  fast  for  direct  inclurion  in  programs  which  iterate  on  other  parameters. 

Pronram  Listings 

The  Fortran  source  listings  of  the  programs  MISTI  and  MIST2  in 
Appendix  C have  been  used  extensively  by  the  author  and  appear  to  be  quite 
robust.  The  author  will  be  most  appreciative  of  comments,  suggestions  and,  in 
particular,  description  of  any  errors  that  might  be  found.  The  first  listings  are  of 
the  two  main  programs  MISTI  and  MiS'n,  followed  by  the  subroutines  In 
alphnbetical  order. 

Proeram.s  Nnl  l.istwl 

Three  of  the  required  subroutines.  DEL2,  DQG32  and  DZERO  are  ftom  the 
open  literature.  These  programs  are  copyrighted,  and  hence  are  referenced,  but  noi 
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listed.  Two  of  the  unlisted  progiams,  DEL2  &nd  DZERO.  requite  itiodiRcelion. 
The  required  changes  are  lisied  in  Appendix  C. 


CHAPTER  6 
CONCLUSIONS 


Accuracy 


Static  Snitinon 

The  mictostrip  equations  derived  in  Ihis  dUsenatiofl,  and  iraplemeiited  in  the 
compuler  progiams,  are  thought  to  be  at,  or  near,  the  accuracy  limit  for  static 

0.001  Z and  almost  entireiy  attributable  to  the  omissian  of  refraction  effects.  The 
derived  equations  are  eiact  for  free  space  and  at  the  limit  for  high  dielectric 
constant.  In  between.  Ihc  effect  of  Uie  series  component  follows  bound  charge 
which  is  praporiional  to  (1  • 1/eJ.  This  is  a derived  relationship  based  on  the 
assumption  that  the  series  field  component  is  uniformly  distributed  above  the 
parallel  component  to  eliminate  refraction.  This  assumption  properly  includes  all  of 
the  eRetls  of  the  series  dielectric,  with  liie  exception  of  variation  of  the  static  field 
distribution  in  the  series  region  for  intermediate  values  of  e,.  The  parallel 
component  of  capadlance  is  exactly  characterized  by  the  equations. 


A further  implication  of  high  accuracy  can  be  seen  from  the  ratio  of  eatemal 
U>  internal  capacitance.  For  alumina,  over  92  per  cent  of  the  field  is  entirely  in 


dielectric 


nbalanced 


strongly  implies  that  propagation  il  accumely  characlerized  by  the  slatio  equalions 
for  alumina.  For  lower  dielectric  coiulant  a greater  proportion  of  the  fleld  is  in  the 
mined  dielectric  (series)  r^on  bur  the  disparity  in  velocity  between  the  intemal  and 
eitemaJ  Fields  is  also  smaller.  These  effecis  are  analyzed  ui  Chapter  5 to  obtain  a 

a wave  solution.  This  analysis  indicates  that  the  maximum  di^eience  from  a wave 
sotudon  should  be  in  Ihe  neighborhood  of «,  = 23.  Unfortunately,  the  analysis 
cannot  provide  absolute  differences  but  these  are  thought  to  be  small. 

Future  Work 

In  years  past,  the  author  has  done  extensive  experimental  work  on 
generalized  microstrip,  particularly  on  applicabons  relating  to  Ihe  design  of  baJun 
transformers  but  using  less  nccurale  design  methods.  There  is  a need  for  updated 
experimental  vehftcabon  of  the  accuracy  esUmates  for  the  theory  of  this 
disseriabon.  particularly  with  reaped  to  wave  propagation.  There  is  also  a need  for 
rtevclopment  of  more  comprehensive  and  efficienl  de&gn  software  for  generalized 

For  some  applications,  c^»citive  coupling  (DC  blockagel  is  required  for  use 
of  the  balun  transformer.  An  investigation  of  coupling  methods  for  tapered  line 
balun  transformers  could  lead  to  improvement  in  performance  for  a variety  of 
applications. 

The  direct  application  of  balun  imnsformers  as  feeds  for  spiral  amennas  and 
for  doubly  balanced  mixers  have  been  the  principal  applications  for  generalized 
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microslrip.  A potential  area  for  further  application  is  to  the  design  of  componems 
in  balanced  microstrip  using  baluns  for  input  and  output.  Components  such  as 
branch  line  couplers  can  be  designed  in  balanced  microstrip,  possibly  with 
considerable  benefiL  As  an  eaample,  uring  a 0.010  inch  alumina  substrate,  the 
conductor  width,  for  100.0  Q unbalanced  microstrip,  is  O.OOI4I  inches  whereas  for 
balanced  microstrip  it  is  0.00496  inches,  or  nearly  5 limes  as  wide.  At  higher 
impedance  the  difference  is  even  more  pronounced.  At  180.0  II  unbalanced 
microstrip  gives  a conductor  widlh  of  0.000065  inches,  virtually  unrealizable, 
whereas  balanced  microstrip,  for  Ihe  same  impedance,  has  a widlh  of  0.00104 
inches,  over  IS  times  wider  than  its  unbalanced  counteipan  and  nearly  as  wide  as 
unbalanced  microstrip  at  100.0  I).  This  raises  the  possibility  of  fabricating 
balanced  branch  line  couplers  with  more  sections,  and  ihus  more  bandwidth,  than 
could  be  achieved  in  unbalanced  microstrip.  An  invesligation  of  design  applications 
for  balanced  miciosirip  could  lead  to  significant  developmenls. 


The  careful  reader  has  probably  noted  that  Ihe  basic  theory,  used  here  for 
geoetalized  microstrip,  can  be  applied  to  other  problems  involving  conducting  strips 
on  a dielectric  sheet.  As  an  example,  coplanar  configuralions  are  of  significant 
interest  and  diould  be  readily  unalyzable  by  the  methods  of  this  dissertation. 


APPENDIX  A 

CONFORMAL  TRANSFORMATION  OF  GENERAUZED  MICROSTRIP 
The  exKl  confonnal  tnuisfomiaiions  of  generalized  mierosirip  result  in 
equations  involving  specialized  inl^rals  and  functions.  For  those  tinfamiliar  with 
conformal  mapping,  the  Erst  section  Is  a brief  outiine  of  confomial  transformations, 
the  Schwarz-CItrisioffel  transformation,  and  the  special  integrals  and  funcbons  that 
arise  in  the  derivation  of  the  transformations. 


Consider  the  relationship  between  the  z and  tn-planes: 


Z = F(w) 


where  z=x  + i>andw=B  +iv.  If  F is  an  analytic  function  it. , satisfies  ihe 
CauchyRienumn  conditions,  3x/du  = 3y/3v,  dx/Sv  = -3y/3a.  then  Ihe 
transformation  is  conformal.  Angles  are  preserved  by  the  transformation  and 
squares  in  the  w-plane  map  into  curvilinear  squares  in  the  z-p!ane.  To  obtain  this 
transformation  the  function  F can  be  chosen  at  random,  provided  it  is  analytic,  or  it 
may  be  found  by  the  Schwaiz-Christoffel  method  which  guarantees  it  to  be  anaiytic. 
Separation  of  z = F(w)  inlo  real  and  imaginary  pans  gives  the  two  equations, 

X = f,«t.v)  and  y = fjfn.v) 


V - constanL  which  form  an  orthogonaJ  map  in  the  z-plane.  In  some  cases  the 
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inverse  transfOrmalion 

« = fidj)  and  V = rj(Jtj), 

can  be  found,  permitting  bi'directiona]  mapping  of  points  between  the  two  planes. 
In  those  cases  where  the  inverse  liansformaticm  cannot  be  found,  bi^liiectional 
mapi»ng  can  usually  be  <hMe  iteratively. 

Conformal  uansfoimations  have  been  widely  used  in  the  study  of 
two-dimensional  field  problems,  particularly  in  electrostatics.  Both  stripline  and 
microstrip  are  generally  well  suited  by  their  inherent  geometry  for  analysis  by  the 
Schwaiz-Chrisutffd  method.  This  is  not  always  the  case,  however,  and  some 
important  microstrip  configutations  cannot  be  directly  transformed  by  the  method. 

The  difficulties  arising  in  the  Integration  of  the  Schwara-Christoffel 
polynomial  for  compliceled  polygons  are  well  known.  The  number  of  slgnificanl 
geometries  that  have  been  solved  in  closed  form  since  around  1S70,  when  Schwarz 
and  Chrisioffel  published,  probably  number  no  more  than  a few  hundred  and 
possibly  much  less.  For  those  cases  where  closed  form  solutions  cannot  be 
obtained,  numerical  integration  or  series  expansion  of  the  transformation  can  be 
used.  These  methods,  hovrever,  lose  many  of  the  advantages  of  the  closed  analytic 
solution  and  there  are  usually  better  numerical  methods  available  than  numerical 
miegration  of  the  Schwarz-Chrisloffel  polynomial-  This  difficulty  in  obtaining  a 
closed  form  analytic  solution  is  one  of  Ihc  primary  drivers  of  research  into  modern 
numerical  methods.  Urgent  problems,  usually,  cannot  wail  for  an  analytic  solution. 


been  formulated  In 


They  provide  not  only  design  equations  but  also 


for  numerical  methods. 

Restrictions 

Conformal  mapping,  including  the  Schwarz-Chiistoffe!  transrormadon,  is 
restricted  to  two  dimensions.  The  potential  solutions  obtained  are  those  for  the 
tsvo-dimensional  laplace  equation  and  are  exact  only  for  purely  transverse  fields 
and  homogeneous  dielectric.  For  purely  transverse  propagadng  waves  f.e.,  pure 
TEM  waves,  the  two-dimensional  Laplace  soludon  is  also  a solution  for  the 
one-dimensional  wave  equation.  Applicadon  of  the  Laplace  soludon  to  waves  for 
the  inhomogeneous  dielectric  case  requires  the  quasi-TEM  assumption  which 
principally  involves  assuming  that  there  are  negligible  higher  order  TE  or  TM 
modes  and  n^ligible  losses  or,  staled  more  generally,  that  there  be  n^ligible 
longitudinal  field  from  any  source.  Refraction  at  the  dielectric  boundaries  cannot 
be  treated  by  conformal  mapping  which  restricts  solutions  to  homogeneous 
dielectric  unless  refracdon  effects  are  small.  In  addidon  to  these  restriedons,  there 
is  a pracdcal  limitadon  on  the  number  of  right  angle  vertices  in  the  polygon.  This 
is  due  to  the  malhemaiicai  difficulty  of  integrating  the  resulting  polynomial 
equadon.  To  minimize  the  number  of  vertices,  conductors  are  usually  assumed  to 
be  thin.  The  following  summarizes  the  usual  assumpdons; 

1.  The  longitudinal  electric  field  is  negligible. 

2.  The  conductors  are  immersed  in  a homogeneous  dielectric. 

3.  Conductor  and  dielectric  losses  are  negligible. 


4.  Conduaor  thickness  is  negligible. 


The  Schwatz-Chrisloffel  imnsfonnation  maps  the  interior  of  an  arisittaty 
closed  polygon  onto  a half-plane.  There  are  certain  formal  restrictions  on  the 
nature  of  the  polygm  but  within  those  restrictions  the  transformation  can  always 
obtained  in  polynomial  form.  Once  the  polynomial  function  is  obtained,  the 
transform  fitncticm  can  be  discovered  by  integration  provided  that  it  can  be 
integrated.  The  procedure  for  obtaining  Ihe  Schwarz-Chrisioffel  derivative  can  be 
staled  in  the  following  rules: 

1.  Traverse  the  polygon  perimeter,  tnainlaining  the  interior  on  the  left. 

2.  Note  the  internal  angles,  a,,  al  each  of  the  n venices. 

3.  The  eiponents,  p,.  of  the  resulbng  polynomial  function  are  given  by 

4.  The  sum  of  the  external  angles,  after  traversing  the  entire  periphery,  miis 
be  2t  or.  equivalently,  the  sum  of  the  exponents  must  be  -2. 

5.  If  the  polygon  is  in  Ihe  a-plane  and  ihe  r-piane  is  the  transform  half-plan 

where  Z)  is  a conslam.  possibly  complex. 

Ellintic  Integrals 

The  natuie  of  the  transform  functions  arising  from  integration  of  the 
Schwaizian  derivative  depend  on  the  angles  of  the  polygon's  venices.  Angles 
■ ±t/2,  right  angles,  result  In  exponents  p.  = f 1/2,  the  square  root  of  the 
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elliptic  uilegrals.  There  are  exceptiona  but  in  general  the  following  rules  will 
correctly  predict  the  nature  of  the  iransfonn  function: 

1.  Up  to  two  right  angles  result  in  elementary  functions. 

2.  Three  or  four  right  angles  result  In  elliptic  integrals. 

3.  Five  or  rrrore  right  angles  result  in  hypetelliptic  integrals. 

involving  elliptic  integrals  and  elliptic  Ivnctions.  For  those  unfamiliar  with  these 
iniegtals  and  lunctions  the  following  descriptions,  notation  and  definitions  wiU  assist 
in  following  the  derivations. 

There  are  three  types  of  elliptic  integrals  designated  as  first,  second  and 
tiiinl  kind.  Additionally  each  can  be  incomplete  or  complete.  Notationaily  these 
are  usually  written  as  shown  below. 

BlipUc  Inicgcil  of  the  Firsi  am 

Ffw.h)  - incomplete  elliptic  integral  of  the  first  kind,  with  amplitude  and 


n?,*)  n /— 

The  first  form  of  the  integral  is  known  as  the  1 
L^endre  form.  The  Legendre  form  is  obtained 
substitution  r — sin  Elliptic  integrals  of  the 
to  be  complete  when  the  upper  limit  of  the  Legt 


(A-1) 


by  the 

lird  kind  are  said 


V “ t/2. 


Ffrtlji)  = K{k)  ■ complete  ellipUc  inlegral  of  Uie  firsl  kind,  with  amplitude  t/2 


and  modulus  k,  usually  written  K. 

Bliotic  Inteeial  of  the  Second  ICind 

E(v,k}  - incomplete  elliptic  integral  of  the  second  kind,  with  amplitude  ^ and 
modulus  k. 

4 

« yi  - {•  « 

E(t/2J()  = E(k)  ■ complete  eiliifflc  int^ral  of  the  second  kind,  with  amplitude  »/2, 
and  modulus  k,  usually  written  E. 


nCv.Q^.t)  - incomplete  elliptic  inlegral  of  the  Third  kind,  with  amplitude  v 
parameter  o’  and  modulus  k. 

n(,..».«n/ ■/ 

n(tr/2,oi',k)  - complete  elliptic  integral  of  the  third  kind,  with  amplitude  t/2, 
parameter  o^.  and  modulus  k. 


The  Modulus 

The  modulus  k can  take  any  real,  imaginary  or  complex  value.  In 
engineering  applications,  however,  it  is  usual  lotakeOskSl.  The 
complemcnisuy  modulus  k"  is  lelaied  to  the  modulus  by 

Complete  elliptic  integrals  of  Uie  firsl  and  second  kind  with  modulus  k'  are  usually 
wriuen  K'  and  F. 


The  upper  inc^ratio 


the  ampliluile 


imesrals.  Like  the  modulus,  it  can  be  real  imaginary  or  complea  but  are  uauaJly 
taken  0<rsland0<srS  t/2.  Notationally,  am  u • ani(u,t)  = [u  real] 
and  am"'(«>,*)  = ani"'(sin"'/,k)  = FX^.k). 


Hk  inverse  of  the  elliptic  integral  of  the  fits!  kind  was  investigated  by 
Jacobi  and  Abei.  Jacobi  caiied  Ibis  inverse  function  liie  "sinus  amplitudimis,"  now 
more  commonly  called  tiie  sine  amplitude  u anti  written  sn(n,k)  or  sn  u.  The 


eUneti  by  the  equation 


th  the  circular  limci 


u = sn-'(r,*)  = /S(sin-‘r,«  - F(vJ0 
and  the  corresponding  direct  relationships  are  defined  by 

I = sn(u,*)  = sin  V,  and  ^ = am(n,k). 

hyperbolic  functions,  an  imaginary  period.  The  real  period  is  AK  and  the 
imaginary  period  is  2iX'.  Two  additional  elliptic  functions  can  be  defined 


(A-4) 


nplitude.  The 


(1  - These  have  periods  (iK,  2K  + 2JT')  and  (2/f,  dUT).  respectively. 

Additional  elliptic  Functions  are  dellned  in  terms  of  the  reciprocals  and  quotients  of 
Ihese  three  basic  functions,  giving  twelve  elliptic  functions  in  ali.  In  the  derivations 

= tnfu.ih)  = sn  u/cn  u.  will  be  used. 

The  elliptic  functions  are  usually  written  sn  u,  co  u.  dn  u and  In  u rather 
lhan  sn(u,it').  cnfu.il),  dn(ii,il),  and  in(u,jh),  If  Ihc  argument  is  v,  ihey  become  so  v, 
cn  V,  dn  V,  and  In  V instead  of  sn(v,jh’),  cn(v.i'),  dn(v,i'),  and  tn(v,i').  Specifying 
the  modulus  in  these  cases  is  unnecessary  since  k always  goes  with  u and  V with  v. 
Notation 

The  notation  used  for  elliptic  integrals  and  functions  is  notoriously  chaotic  and 
unstandardized.  For  this  reason  a single  reference  was  selected.  'Byrd  and 
Priednian  [15].'  At  appropriue  paints  in  the  following  text  the  notation  ‘B.F. 
niui.m'  appears  following  an  equation,  meaning  that  the  source  of  the  equation  or 
identity  is  Byrd  and  Friedman's  equation  Departures  from  their  notation 

will  be  noted  where  Ihey  occur.  An  alternative  reference,  which  may  be  more 
generally  available  lhan  Byrd  and  Friedman,  is  the  'Table  of  Integrals.  Series  and 
Products'  by  Gradshieyn  and  Ryzhik  [16].  Their  treatment  of  elliptic  inlegials  and 
functions  is  quite  thorough,  though  not  as  complete  as  Byrd  and  Friedman. 


The  z-plane  lo  w-plane  transfonuUon  is  obBincd  through  the  use  of  i 


imennediale  r-plane  Uansrormalion  with  branch  points  common  to  both  the  z and 
w-planes.  This  tnnsformadon  was  obtained  by  Black  and  Higgins  |2]  using  values 
at  die  branch  points  different  from  those  used  here.  Their  solution  can  be 
transfonned  to  obtain  the  result  derived  here;  It  is  useful  and  illustrative,  however, 
to  carry  out  the  entire  derivation  independent  of  their  prior  solution. 

The  three  iransform  planes  are  illustrated  in  Fig.  A-l  and  the  transform 
points  are  tabulated  in  Table  A-l.  Note  that  the  sum  of  the  eiponents  is  -2  for 

In  the  process  of  deriving  the  transforms  of  generalized  microstrip,  a 
number  of  definitions  and  identities  relating  to  elliptic  integrals  and  functions  will 
be  required.  All  of  these  can  be  found  in  the  suggested  references;  For 
convenience,  however,  the  applicable  relationships  have  been  collected  in  Table  A- 
2.  These  formulas  are  numbered  l■aa  and  will  be  referred  lo  by  number  as  needed 

The  w to  r-olane  TransformaHon 

Uung  the  internal  angles  of  the  polygon  In  the  w-plane  and  traversing 
counter-clockwise  from  (2)  lo  (7)  and  back  lo  (2).  in  Fig.  A-l,  gives  the 
Schwardan  dtfferenbal. 

(2)  (4)  (5)  (7) 

ehv  = D‘(i  + If*)-''*  <f  + I)-'**  (I  - 1)-!'*  (r  - 1/*)-''** 


Fif.  A-l  Schvan-Chrlstof/el  truu/ormation  diagram  lor  i,  U and  v-planaa. 


Tabulation  of  tranalorm  polnU. 


Point 

® 

@ 

@ 

® 

® 

@ 

Z'plano 

Ini. 

-h 

-h+ib 

-h 

tbi-la 

+h 

t-plana 

-1/Slt 

-lA 

-C2/2 

-1 

♦1 

*Cl/k 

+l/k 

Szponent 

-8 

-1/2 

♦1 

-1/3 

-1/3 

+ 1 

-1/2 

UgtliT 

-E+ir 

-K+lb' 

-K 

rK 

eEtla' 

.»E+ir 

t-plana 

-1/Rk 

-i/it 

-C2/k 

-1 

♦ 1 

♦Cl/k 

♦ lA 

Exponent 

-1/2 

-1/3 

-1/2 

-1/3 

VhS  fSTcSSSs 


= 0'((^  - l/i’)"*'’  (1^  - I)*’'’  * 
- D’/t  -1)''"  {/“  - !)•''’*. 


= D(I  - 

•“’fe 


where  D and  Nirt  consonU.  possibly  complex,  which 
posilion  of  the  iwpLuie  polygon.  Equation  (A-3)  is  set 
(A-I)  with  solulion 


At  point  (4)  r - -I,  w = -m  and  ^ = sb-'(-l)  giving 
w = DF(-rl2.k)  + JVor 
-K(k)  = -DK(k)  + N. 

At  point  <5)  f = +1,  w = +AT(i)  and  ^ = sin-'(+l)  giving 
w = DFl+T/2,k)  + Wot 


(A-6) 


(A-7) 


+K(k)  = +DK{k)  + N.  (A-8) 

Addition  of  (A-7)  and  (A-8)  gives  W =0  and  subtraction  D = \.  With  these  values 
(A-6)  becomes 


w “ fXsin"'/,*). 


(A-9) 


Using  (A-9)  and  (A-4), 
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C»a»ciBiict  Per  Unil  Lani>-  ■",< 

The  free-space  capacitance  per  unit-length  of  generalized  mlcrostrip  can  now 
be  obtained  directly  from  the  w-plane  map  aa  permittivity  times  conducicpr 
widih/separaiion,  Cg/2  = Farads/meter  for  the  upper  half-plane,  or 


Farads/meter  for  the  total  capacitance,  including  the  lower  half-plane.  The 
propagation  velocity  on  the  mlcrostrip  is  v,  = ((«)“’-  For  free  space  v,  - 
(figCo)''*  ~ c,  the  free-space  velocity  of  light.  The  characteristic  impedance  of 
microstrip  is  defined  by  Z = (L/C)"^.  For  free-space  I,  = fig  and  C = Cq,  giving 
2,  - WCo)”’  «•  „ 


lb  . ilS  . J76.7  Q . 


as  far.  eipressions  for  the  free-space  ca 
generalized  micToslrip  have  been  obtained.  It  rer 


(A-12) 

.nee  and  impedance  of 
to  determine  the  relationship 
lucting  srips  by  finding  the 


The  f 
This  polygon 
oblairted  by ' 


transformation. 

The  w to  r-plane  Transformation 
nite  portion  of  the  z-plane  polygon  is  shown  at  i 
includes  infinity  as  verm  (I).  The  Schwarzian 
raversing  the  periphery  of  the  polygon  from  <1) 


etop  of  Fig.  A-1. 
differemial  is 


(2).  (3),  (4).  (5).  (6)  and  (7)  to  (1)  ■ 


the  left,  lunine  through  an  angle  of  2t  radians  between  <1)  right  and  (I)  left. 


transform  points  are  given  in  Table  A-l  and  result  in  Schwaraian  derivative, 


The  derivative  of  the  tv  to  r-plane  transform  was  derived  in  the  previous  section,  it 


Dividing  (A-14)  by  (A-13)  yields 

ikidw  = (lUldtWvIdi)  = D(r  - C//i)(r  + C2/i)(r  -H/itt)-*  (A-IS) 
where  D = D7D‘.  The  left  side  of  this  (A-IS)  is  a function  of  w whereas  the  right 
»de  is  a function  of  r but  (A‘14)  was  previously  integrated,  resulting  in  the 
transformation  r = sn  w.  Substitution  of  sn  w for  t.  the  dummy  variable  ( for  w, 
and  integrating  both  ddes  of  (A-15),  gives  the  integral  that  must  be  evaluated. 


dzidi  = D'(r  - Cllk)(i  + C2/*>(r  +l/flk)-’ 
x((r+  l/i)(r  - !/«(/  + l)(r  - 1))-''*. 


<A-13) 


dw/di  = D'((r  + !/*)(;  - l/*)(r  + l)(r  - I))"''*,  (A-14) 


(so  £♦!/«)= 


(A-16) 


Expanding  the  numerator,  (A-16)  can  be  written 


(A-17) 
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Direct  [niegratioti 

The  5itt  auempt  is  dired  integration  of  (A-17)  term  by  term.  The  integraJ 
of  the  third  term  of  (A-17)  is  found  in  Byrd  and  Friedman  (B.F.  330.51,  330.50, 
and  361.58)  which  are  shown  below. 

B.F.330.51  •’-I’fz — i V; 

{of  - l)(oft^  - 1) 

« n(<p,l/o?yl)  - mji) » (1  - ,« 

^ cnadnu  j 

o,  j,(o, -r  gn  u) 


^.2(e^-l)(a»-t^t(l-gW»)(«^*ttU>-2h») 


2l((a>-IX«’-ik‘) 


This  is  the  hyperbolic  case  of  the  elliptic  integral  of  the  third  kind,  since 
0 S o’  S *“.  The  compienity  of  this  integral  is  such  that  it  appears  unlikely  that 
the  drat  two  terms  can  be  diiecily  integrated. 


It  is  noted  that  the  transformation  must  be  of  the  form 


tfni)  = 0(w)  [ ,.j  + ♦(»,«) 


(A-18) 


i^ined  funclic 


of«i 


ew  it 


when  6(hO  and  f'iR.y)  are  iindelerm 
fuiKlion  or  w that  remains  when  ff  - 0.  With  J!  = 0,  the  tnuisrormalion  reduces 
to  that  for  equal  strips  (balanced  microstrip).  Taking  R as  zero  puts  the  pole  in  the 
r-plane  (l/JU)  at  Infinity,  thus  removing  It  from  ihe  integral.  The  resulting  reduced 
integral,  widi  Cl  ^ C2  = CO.  is 

iMIi-o-O'  I j + «.  (A-19) 

Ihis  equation  Is  mtegrated  using  B.F.  310.02,  yielding 

«w>  [ ,.o  = D'((l/*’)(w  - £(H9)  - (ClVi)V)]  + AT, 

- 0[(l  - af)w  - £(w)]  + IV  (A-20) 

with  D = 

llie  constants  D.  N.  and  02  are  evaluated  by  setting  r and  sv  to  appropriate  values 
at  poinu  (2)  throu^  (7)  using  Fig.  A-1  and  Table  A-I. 

At  point  (5)  z = +A,  w = +K  giving 

+h  » 0(+(l  - af)K  - E(K)]  + N.  (A-2I) 

At  point  (4)  z = -h.  w = -K  giving 

-h  = £)[-(!  - CfilK  + £(A53  4 M (A-22) 

Addition  of  (A-21)  and  (A-22)  gives 

0 = 0 4 2^,  or  iV  = 0. 

At  point  (7)  z = 4A,  tv  » a:  4 iA('  and  using  1-24  to  expand  &K  4 iA'),  gives, 

4A  = 01(1  - Ctf)K  4i(l  - OT’j/f  - £ - i(Al'  - £•)].  {A 


\-23) 
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-A  = 0[-(l  - CtfiK  4i(l  - a^K  *■  E-  HK-  - £’)J.  (A-24) 

Adding  (A-23)  and  (A-24)  gives 

0 = i 2D[(1  - a^)K‘  - K'  - £•], 

(1  - aj*)r 


Subs&niiion  of  (A-2S)  imo  (A-20),  with  N = 0.  gives 


(A-2S) 


(A-26) 


At  point  (5)  t = +h,  w = +AT  giving 

(A-27) 


At  point  (4)  z = -A.  w = -K  giving 
Subtraction  of  (A-28)  from  (A-27)  gives 


Using  1-26,  Legendre's  relationship,  the  denominator  of  (A-29)  is  seen  to  be  equal 
to  —w/2  giving 


Substitution  of  (A-30)  into  (A-26)  gives  the 


z-plane  transformation  for  the 


qxcisl  case. 


(A-31) 


The  earlier  attempt  at  direct  integiaiion  of  (A-17).  term  by  leiin,  gave  such 
a complicated  lesull  Ihai  the  approach  was  abandoned.  It  is  unusual  for  physical  or 
geometric  transfonns  to  yield  such  complicated  lesults.  so  it  is  posable  that  there  is 
cancellation  between  terms.  With  this  possibility  in  mind,  assume  that  the  Integral 
has  the  form  of  (A-lg).  Substituting  (A-31)  in  (A-13)  gives 


where  ~t(w,lt)  is  a function  to  be  determined.  The  negadve  sign  on  ♦(w.S)  is  from 
eq.  (A-30).  Returning  to  (A-17),  it  is  noted  lhat  the  denominator  appears  as  the 
square  of  a funcbon.  Thus,  ifi’/ifw  has  the  form  of  the  derivative  of  a quotient, 
given  by  the  well  known  formula 


Taking  U = «w>  and  V - sn  w 4-  lIRk  gives 
T (wj?)  - 

snw4I/ftt 

with  <Hw)  r^lacing  ~t(R,  w)  as  the  unknown  function  lo  be  determined. 
SubsUtution  of  (A-33)  into  (A-32)  gives 


Comparing  (A-34)  to  B.F.  330.SI,  obtained  earlier  in  the  aiienipt  to  direct 
integrate  (A-17),  shows  that  the  last  term  of  the  iniegiai  has  the  same  form 
(A-33)  and  that  -cn  w dn  w is  a likely  candidate  for  ^w).  Stibstitulion  ol 
-cn  w dn  w for  iS(w)  gives  the  w to  j-plane  transformation  as 


ded  that  Hv)  ht 
actly  establishet 


been  correctly  selected.  The  t 


necessaiyr  to  take  the  derivative  of  (A-35).  equate  it  to  (A-17)  and  delenn 
eonstanu  Cl  and  C2  in  terms  of  the  coefficients  of  (A-35).  The  derivativ 
(A-35)  is  obtained  using  B.F.  314.02,  731,01,  731.02  and  731.03,  whic 


(A-36) 


To  simpliiy  noladon  and  manipulation 
{ = Vm.  Malting  the  indicated  subst 


Taking  the  derivative  of  both  sides  of  (A-17)  gives 

A - p[  (C2/4  - Cllkis  - ClC2lk>] 


Equaling  (A-37)  and  (A-38)  al 


= D[^  + iC2lt  - Cllk)s  - ClCZHf). 
To  establish  the  equality  (A-39)  Ihe  two  sides  must  be  in 


(A-39) 


is  done  by  manipulating  the  left  hand  sde  to  obtain 

LHS  = Hcoeff,  X sn’  w + coq^  x sn  w + cwj^  ].  (A-40) 

Obtaining  the  form,  (A-40)  is  a necessary  condition  for  (A-3S)  to  be  a solution.  In 
the  following  steps  Ihe  LHS  is  manipulated  using  identides  1-6  through  I-IO  from 
Table  A-Z. 

LHS  - U(s  + £)’M  + (s  + 6>’d’  - - s*d^  - - {sd*  - cV) 

= L[(s  + {)=*/  + (s  + - *=s’c*  - (o’  + s*)d*  - {si’c'  - {sd’l 

= L[(s  + (fU  + s»d*  + 2{sd’  + {V  - i^sV  - d’  - {stV  - {sd’l 

= i[(s  + {)’«  4 $*6’  - AV  4 {sd’  - {st’c’  4 e’d’  - d’j 

- t[(s  4 {)’«  4 fV  + i”(s  4 {’d’  - d’l 

= L|(s  4 {)’«  4 (1  - A’ls’  4 (1  - *’){s  4 ({’  - Dd’] 

- LI(S  4 i)’M  4 s’  - {’s’  4 (1  - {’){s  4 ({’  - i)(i  - i’s’)] 

= £((s  4 {)’M  + ,t  - J4,s  + (1  - *S|{j  + {’  _ t>5»sS  _ I + *5st| 

= L(s’M  42Af£s  4 M£’  4 s’  4 (I  - i’){s  4 {’  - H‘(V  - I] 

= i[(l  4 M - *’£’)s’  4 (1  4 2M  - {’)£s  4 (1  4 A0{’  - IJ  (A-4I) 

SubsUludon  of  (A-41)  for  the  LHS  of  (A-39)  gives 

t[(l  4 « - {’{V  4 (1  4 ZjH  - *’){s  4 (I  4 Af){’  - 1) 

= D[s’  4 (C2/i  - CM)i  - C/C2/if).  (A-42) 

The  equality  (A-42)  holds  if 


t(l  4 M - ie{’), 


(A-43) 
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D(C2  - ayk  = Ml  + 2M  - **)£  ana  (A-44) 

- DCianf  = i«’  + «{’  - 1)  <A-*5) 

which  can  be  equivalently  written 

i.  = D/(l  + M - *’{’), 

(C2  - Cl)  = <1  + 2M  - i’ltf/Cl  + M - and  (A-46) 
CiC2  = (!-{’-  Mj’)i’/{1  + M - t’j’).  (A-47) 


Substitudon  of  M = B/K  - 1 and  ( = 1/At  into  (A-46).  and  (A-47),  and  after 


C2.c,~BMLzK^£E1,  (A-48) 

A’£'  - Jf' 

CJ«C3“  _ (A^9) 

R'B  - A’ 

Back  substitution  of  (A-43),  (A-44)  and  (A-45)  into  Ihe  LHS  of  (A-42)  esiablishes 


(A-50) 


Evaluation  of  the  constant  N completes  the  derivation  of  the  transform  i 
The  constant  N is  delermined  by  evaluation  of  (A-50)  at  corresponding 
i and  w-planes,  using  Fig.  A- 1,  Table  A-1. 

At  point  (5)  a - -fS  -HO  and  w = Agiving 


. 2A£Lr  _ 


At  point  (4) . 


It  [ r ' m-j»r+i/«J 

Using  identities  M2  to  14,  1-17  to  19  and  1-29  in  (A-51)  and  (A-S2)  give 


- * - I - ( - 1 ) AJ(«  - £(*)  j -f  W . (A-54) 

Adding  (A-53)  and  (A-S4)  gives  0 = 2iV  or  \ = 0.  This  completes  the  derivation 
of  the  IV  to  z-tiansfbrmalion  which  can  now  be  written  in  final  form  as 


Etguatian  (A-S5)  Is  in  comples  form  with  : = x + iy  and  w = u -t-  iv,  a 
form  dial  is  not  computallonally  useful.  It  can  be  reduced  to  real  form  by 

teal  pan  of  Ihe  right  hand  side  of  (A-SS)  and  y to  the  imaginary  pan.  Direct 
application  of  this  procedure  results  in  a rather  complicated  and  cumbersome 

of  (A-IS).  With  this  procedure  Ihe  equations  for  *<o,v)|».,  = 91[e(u.v)]  and 
y(u,v)  Id.Q  — 5 [0(u,v)]  where  31  and  0 designate  the  real  and  imaginary  pans 
respectively.  The  general  Iransfonnations  are  Ihen  obtained  by  adding  9?(*(u,v,«)| 
and  Uf-KUiV./?))  to  the  result  obtained  for  S = 0.  This  gives  xfu.v.If)  = 9I(e(u,y)] 


andXir,v,B)  = a|0(u,v)]  + a('f-(u,v,«)l. 


(A-56> 


-)  l„  - - SCv.*- ) + j . (A-S7) 


. cdd,(c,»  - etWh  B = «,c,(*V  + (I’d,*),  C - sd„  D = Cd5,c„ 
= 1 - i\‘.  wilh  s = sn(u,k).  c - CT(u,k),  i = dn(u,li),  s,  = sn<v.k'). 


X = xlh  Hid  r = ylh.  I 
special  case  J?  = 0.  gi 


(A-63) 
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The  general  transromiabon  Is  given  by 


(A-64) 


« ’ A + J' 


Inverse  Transfomalions 

The  possible  exisience  of  inverse  transformaiions  of  the  form,  ii  = fiCJ,y) 
*!'<!>’  = Jitro').  briefly  discussed  in  section  Al.O,  Unfonunalely,  for 
generalized  microstrip,  no  such  inverse  exists  and  inversion  must  be  carried  out 
iteratively.  For  the  special  case  11  = 0,  inverse  transformations  can  be  obtained  in 
useful  form,  though  not  in  the  ideal  form  ^own  above.  The  derivable 
transformations  are  v = and  n = u(y,v,*').  The  generality  of  these 

transformations  is  somewhat  limited  by  the  appearance  of  a as  an  independent 
variable  of  v and  v as  an  independent  variable  of  u.  Decile  this  lack  of  generality, 
they  do  provide  a means  of  directly  mapping  the  dielectric  boundaries  from  the 
a-plane  to  the  w-plane.  The  derivation  of  v = v(».o,ifc)  begins  with  equation  (A- 
62).  which  can  be  written 


(A-66) 


* 5**'  *“'’**'•“'*  /’into  (A-56)  and  solve  fbr  s,^, 


Pd'  * *’««/ 

. iW-«V 


v{X.B.i)-F(r.,i'). 


i r.-s 


(A-69) 

(A-70) 

(A-71) 

(A-72) 

<A-73) 


The  iiivene  iransrornialion  u - u(y,v,l:‘)  is  oblniited  from  equation  (A-63) 
by  solving  for  d'.  giving 

(A-74) 


M-ger 


^ Qdj  - ih'a,ejdj 


{A-75) 

(A-76) 


(A-T7) 


' the  expression  for  alh  is  obtained  by  evaluating  (A-S5) 


point  (6),  where  from  T^le  A-1,  z = A + ij,  w = *:  + w’and/  = CUk. 
Substitution  of  these  values  for  z,  and  w in  (A-SS)  gives 


cn(Jir*ifl')<lii(Arvifl’)| 


h-78) 

1-79) 


From  (A-10),  r = sn  w;  at  point  (6),  f =■  C//i  and  w = AT  7-  ia'  giving 

Ci/li-sn(A:*i<»')  (A-80) 

Using  formula  1-21,  vrith  u = ff,  v = <r‘,  and  aRer  some  manipulation,  it  is  found 
that  sn(/f  -I-  ia')  = l/dnla',*")  giving  CUk  = l/dn(o',i')  or 

(A-81) 

Identity  1-8  gives  an  expression  for  sn(a',b')  in  terms  of  in(a',k')  which  can  be 


(A-S2) 


larly,  using  identily  1-9,  cn<a',*')  ca 


cn(a'.f)-j 


1-83) 


Reduction  of  cn(K  + id')  and  dn(K  + id')  is  done  through  use  of  identities  1-23 


and  1-23,  giving 


ca()f*i<i')--i| 


[fao(d'.if) 

[ dn(a',f) 


(A-84) 

(A-85) 


Subsdtudon  of  (A-Sl)  and  (A-82)  into  the  right  hand  side  of  (A-S4)  ^ves 

cn()r*io')--l  - ****** 

SimUarly  subsdtudon  of  (A-81)  and  (A-83)  into  the  right  hand  side  of  (A-3S)  gives 
do(g:*ifl')-i/l-C/^  **•*’> 

Subsdtudon  of  (A-80),  (A-86)  and  (A-87)  inio  the  last  term  of  (A-7S)  gives 

gtfgtiii'jdnfgt-ifl’)  _ , |i/CJ^  - - C/^|  (A-88) 

saiK*io')*VRk  [ CI*]/R  J' 

The  Imaginary  part  of  die  first  two  terms  of  (A-78)  Is  given  by 

a|(|l  - IXAT-r  id')j . i (,£  - l)fl'  . S£(*:  • la') . (A-89) 


The  imaginary  part  of  E(£  -h  in')  is  evaJuated  using  idendty  1-27  which  gives 

SE(K*  id')  - d'  - Efa'.k’)  * • (A-®0) 

From  I-I8,  cn  £ = 0,  giving 

3£(Ar*io’)  ■ o’  - £(0',*') . (A-91) 

Subsdtudon  of  (A-91)  into  the  right  hand  side  of  (A-89)  gives 


3[(-^-lXArs-ld’)-£(£*ia')]-i[(|^-l)a'-ro'-£(o-,i’)] 


(A-92) 


Substitution  of  (A-92)  and  the  riglil  hand  side  of  (A-88)  into  (A-78)  gives  Ihe 
espiesaon  for  alh  as 

Mil  I (A-93) 

It  I Jf  ^ \J  C1*\IR  \ 

In  (A-93),  am  o'  = ^ is  the  amplitude  and  k'  is  the  modulus  of  the  elliptic 
integrals.  Accordingly,  (A-93)  can  be  rewritten  in  Legendre  formal  using  M,  svilh 
V - a'  = ft*,*'),  giving 

" I g . «,.f ) - I , tA-Xl 

It  [r  ci  + i/R  j 

y arcsin[8ii(a*,jt*)l  ■ arcsip^^^^  (A-95) 


Uang  equations  (A-94)  and  (A-9S)  alh  can  be  evaluated  The  equations,  however, 
are  not  numerically  stable  for  ^ very  near  t/2.  A better  evaluation  of  ihe  elliptic 
integrals  is  obtained  using  the  Bulirsch  algorithms  |17].  The  upper  limit  of 
integration  for  these  algorithms  is  T,  = tan  «>  = tn(a',h')  = sn(fl',*')/cn(a',jt'). 
The  tnfo'.i’)  can  be  written  using  (A-82)  and  (A-83)  as 


r.“<£t*li. 

' ) 

and  o'  = FU,.k-). 
wiUi  \KI2  S a'  S AT). 

10  (A-94),  and  with  some  manipulation,  j 


(A-96) 

(A-97) 


o/A  ■ — |r r,  ,f ) - «r,  ,i' ) ♦ r.i(  i - ci»)  -- 


\Cf 


.(A-98) 


(A-99) 
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By  a similar  procedure  the  eapression  fc 


and  *'  = 

th  1^/2  &b‘  S IT]  and  [b‘  2 a’]. 


Exaression  for  a/b  - b/h  vilh  R = 0.0 

For  the  special  case  of  balanced  microstrip,  R - 0.0,  the  formulas  s 


(A-102) 

(A-IOJ) 


alb\^‘blb-^\^n.T,.k')-EiT,.k'y-< 


ZkililE!!].  (A-U 


The  final  paramelcr  to  be  evaluated  Is  Ug,  |-/f  £ Ug  s 0],  At  point  (1), 


and  w = Ug  + iAT  with  r = sti  IV  giving 

- I/Bt  = sn(ug  + iT).  (A-106) 

The  right  hand  side  of  (A-IDS)  is  evaluated  using  1-21,  giving 

sn(U(,  + iA:’)  = I/(tsn<B|„i).  (A-107) 


Substinilion  or(A-l06)  into  (A-IQS)  gives 

- 1/Jlt  = or 

- ft  = sn(i(o,*)  (A-108) 

Squaring  boUi  sides  of  (A-107)  gives  sn’(n,,*)  = JH,  so  cn*(U(.i)  - I - ft*  and 
tn‘<uo,<b)  - sn*(uo»/cn^iio^),  vriiich  can  be  wriiicn 

(A-109) 


with  -ft  £ u,  £ 0. 

This  completes  the  derivation  of  equations  for  the  free-space  parameters  of 


APPENDIX  B 

ANALYSIS  OF  THE  DIELECTRIC 
POLARIZATION  CAPACITANCE 


The  anfliysis  of  ih 
elearificacion  of  Uie  inne 
Figure  B-l  is  a di 
space  coordinates  (z-plan 
two  half  planes  using  the 


t appendix  is  exact  for  the  dielectric  polarization  due  to 
faces  of  Ihe  conducting  strips, 
gram  of  the  transformations.  The  polygons  representing 
) and  flux-poicnrial  coordinates  <w-plane)  are  mapped  to 
well  known  tiansfbnnation  r — sn(z'^  and  r*  sn{w^). 


Those  unfamiliar  with  this  transformation  may  refer  to  the  section  A-2, 1 of 
Appendix  A where  the  transformation  is  derived.  These  polygons  are  of 
dimensions  N by  N"  uid  P by  P where  N,  ff , P utd  P'  aic  complete  elliptic 
integrals  of  the  first  kind  with  moduli  k and  p,  respectively.  The  z and  z’ -planes 
are  connected  by  the  scale  factor  m = NIn  = If/n’. 

The  inlermediale  half-planes  ara  connected  by  the  bilinear  ttaiisformalion 


At  point  (5).  P = I/p.  I 
(B-l)  gives 


f,-  (/-<fHr-c) 

(d-cX/-r)' 

c = 0 and/=  o».  Substitutii 


(B-l) 


101 


n'eab 

Dielectric  I 

-planoi  (UicrcsU 

mh+iN' 

d) 

1 

Dielectric 

@ 

me  z’-dzU 

N(k) 

t=Sn(z'.k).  Vi 

=N/n=N'/n' 

® 

@ 
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‘■{d-c)(f-t) 
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Dielectric 

Fig.  B-1  TranaformAlion  dlAgram  for  z',  t,  V. 
detormisoUon  of  the  polerization  cepacltuiee 
mlcrootrip. 


and  v-planeo  for 
of  generalized 


£® " -e;® 


W)  » l/[isn(mi>], 


unhM.  yielding 


l/=l-p'»in(B5),  giving 


where  A = toMH.  and  fl  ==  t6/4A. 

The  geometric  capacitance,  computed  as  dielectric  constant  times  conductor 
width  divided  by  conductor  separation  in  the  w-plane  of  Fig.  B-l,  is 


Equation  (B-7)  is  the  total  internal  capacitance  of  the  microslrip  in  the 
z-plane  of  Fig,  B-1.  It  includes  not  only  the  capacitance  due  to  polarization  of  the 
dielectric  but  also  the  capacitance  due  to  free  charges  on  the  conductors.  Because 
of  this,  direct  use  of  (B-7|  as  the  polarization  capacitance  would  overstate  the 
capacitance  by  the  fraction  of  capacitance  attributable  to  flee  charge. 

Dielectric  material  between  the  conductors  of  microstrip  increases  the 
capacitance  by  polarizing  charges  at  the  conductor  faces  that  would  otherwise 
contribute  to  the  internal  electric  field.  The  phenomena  of  polarization  can  be 
visualized  as  the  binding  together  of  dipoles  to  form  a dipole  chain  under  the 
influence  of  the  applied  potential  with  their  free  ends  terminating  on  the  cwiductor 
surfaces.  These  dipole  chains  constitute  the  polarization  held  and  do  not  contribute 
to  the  electric  held.  The  voltage  across  a microstrip  transmission  line  can  be 


(B-8) 


is  the  geometric  capacilai 


total  charge  Q contribute 


tates  that  only  1/e,  of  ilv 
le  conductors.  This 


(B-IO) 

le  bound  charge,  is  neutralized  by  the  di^tric  and  therefore  contributes 
to  the  electric  field  within  the  microstrip.  Capacitance  is  defined  as  C ™ 
using  (B-IO)  and  (B-7)  it  is  seen  that  the  polarization  capacitance  is 


(1--5. 


c - rn 

..^r«  . 


(B-Il) 


APPENDIX  C 

PROGRAM  SOURCE  CODE  USTINGS 


PROGRAM  HlSn.FOR 

PURPOSE 

Che  exeec  fenulaclOTt  of  eq.(3*23).  Clven  Che  fcee-spece 

e Uieleoccic  sheec  (E),  Che  epeclal  pereaecere 


.IinmiV,DFARAH.OQC32. 


C SUBROUTIHZ  AND  FWCTIW  rKOCRAHS  KEIJUIftBS 


C R£FE££HCE 


FASAMnTR(CC-2. 99792AStDO) 
10  HB1TE(S,1; 


CALL  DELO(ZO] 
IF(a:*XK.LT.l.D-40)  THEN 
VE1TE(S,4) 


CALL  DCONST(R,C2,C2) 


CALL  SPASAH(R,C1,C2,AP,BP,UC,AK,EH) 

CALL  SCF(AH,BH,CP) 

CALL  SAREA(R.AP,UO,CPHAX] 

IF(R£PS.LT.l.DO)  GOTO  10 
C1-CP<(REPS • 1 . DO)«AP 

CE-(l.DO-AP)/(l.BO-(l.DO-l.lW/EEPS)»(CPllAX-CP)/<l.I>0-AP)) 

EFFK-CI«CE 

VP.1,D0/DSQRT(EPTK> 

CAP-1. D4/(CC»VP»Z) 

HR1TE(5,12) 

HR1TE(5,7) 

CI-CPPCEEFS-l.DOtAP 

CE-(1.DO-AP)/(1.DO-0.5}6DO»(1.DO-1.I»/REPS)*(CP-AP)/<1.I)0-AP)) 

VP-l.DO/DSqRT(EPnC) 

CAP-1 . D6/<CC»VP*Z1) 

WRITE{5.7} 

WRITE{5.B)  Zl.AH.BK.VP.CAP 
WRITE{5,11)  100.B0»(2-Z1)/Z 
WRITE(S.U) 

C FORMAT  STATEMEHTS 

1 PDRMATC//'  Type  113.5.LE.ZO.LE. B40.0) , (0.0)  Co  quit.  20-  ',$> 

2 FORHATC  Type  (O.O.LE.R.LE.1.0) , l-l.O]  Co  change  20.  R - ',$) 

3 PORHAT(D2S.16) 

A FORHAT(UX'IHFEOANCE  TOO  LOU,  REENTER  20’] 

3 PORHATdBX’IHFEOANCE  TOO  HICK,  BE  ENTER  20') 

B PORHATC  Type  EPS.  [EPB.CE.  1.0| , |0.0|  for  nev  Z.  EPS  - ’,$) 
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7 roiWAT(6X,  ’Z’  ,lix,  ’«/h’  ,13X,'b/h'  ,13X,’Vp'  ,llX,'c«p  pf/B’) 

i roEUUT(SG16.t) 

7 roUUT(/33X'*PFROXIIUTiaN’  > 

11  Ft3a«XT(/10X'  BIFF  » - lEK).*(EXACT-APFilOX>/EXACT  - ’.FIB.IS) 

12  FOUUT(/30X,' EXACT  FOKHUIATION’ ) 

13  F0BU1AT(/) 
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C REHAKKS 


C abioluta  »e«irMy  of  .ppcoilMCly  O.OC11»2,  including 
C tcfroctlon.  A ucrnlng  notice  1.  given  -hen  cuctlde  the  region. 


C SUBKOUTINE  AND  FUNCTION  FROCFAHS  REQUIRES 


C REFERENCE 


CHARACTERel  UE 
INTEGER  lER 

EXTERNAL  D»2_l 


Ill 


IF(ll£PS.t.I.l.DO) 


EFFL  - 0-5DO*(BEPS  + l.DO) 

ML  - DSlJRT(EFn.)»TASaer*0.9MD0 
20U  - DSgRT(EFFH)*TA!IGST*1.001DO 


CALL  DELO(ZO) 
TEST  - KKP/XKKP 

URITE(3,6) 


CALL  DBLO(M) 
«Rm(5,7) 


ZOO  - DZEBO(ZOL,ZOU,DF»2_1,TOL 


112 


WRITE  (S,8) 


WRITE  (S.9) 

WRITE  <5,10> 

W8ITE(5,11)  Z/SU,BU*AH,VP,BU*CAr 
WRITE{5,12)  MO/M 


2 MRtUTtAl] 

3 It}RIUT 

&(//’  Type  Upedsnc*,  [0.0|  ts  quit,  (-l.O]  Co  change  U/B.  Z - ',$> 

fcC  Type  ( EPS  .CE.  l.O  1,  tO.O)  for  new  Z.  EPS  - 
3 F0R1UT<D23.1S) 

6 FORl«T(17X,'e«  2 TOO  HIGH,  REBSTER  Z ♦*»') 

I P0R11\T(I7Z,'»**  Z TOO  low,  RBEtrTER  Z W) 

( PDIUUT(27Z,‘lI)iBALANCED’) 

9  PDIIHAKZSX.'BALAHCED’) 

10  F0RMAT(IX,'Z',UZ,'eAM3X,’vp',llX,'eep  pf/o' ) 

11  P0RHAT(4G1E.6) 

12  FORMAT(21X.’  20  - ^GZO.IS) 


SUBROUTINE  MRE\(R,AP,UO,CPHAX) 

C PURPOSE 

C Incegracee  Che  dielectric  area.  In  Che  w-plana,  of  generaltsed 

C elcroicrlp. 

C CALL  MREB(R.AP,UO,CPHAX) 

C DESCRIPTION  OP  PARAMETERS 

C CAMAR  - RESUII  VALUE  OP  DIELECIRIC  AREA  IN  DOUBLE  PRECISION 
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C RBHAIUU 


• INPUT  VALUE  OP  THE  PARAHETES  E IN  DOUELE  PRECISION 

• INPUT  VALUE  OF  A"  IN  DOUBLE  PRECISICW 

• INPUT  VALUE  OP  UO  IN  DOUBLE  PRECISION 


C SUBROUTINE  AND  FUNCTION  PSOCRAHS  REQUIRED 
C PREVIOUS  SUBROUTINE  AND  FUNCTION  PROCRAH  CAl 


C ElnEinB  rouclnA  DZERO. 
C REFERENCE 


IHPLICIT  REAL'S  (A-K.O-Z) 
EETERHAL  DXOFV,  DFINDU,  DFINDV 


CALL  DSa>(SU.CU.DU,UO*XXX,lKP> 
X-SU/CU 

EU-DEU  CX , XKP , 1 . DO,  XKP'XKP) 


VO  - BZZRO<VL,VH,I«OFV,TOl,ISR) 
VO-O.99W99$D0 


CAIX  OqC3Z<in..UK,BFiniV,AREAl) 


IF(ft.EQ.l.BO>  GOTO  30 


VH-O. 99999 999 DO 


CALL  DQG32(UL,1JS,DFIMDV,AREA5) 

C Skip  Incpgracion  of  AREA3  If  IER-1,  l.o.  no  Inflocclon 
IF(IER.Eq.l)  GOTO  30 


Ul^-0.9999999D0 
UM-UO  t O.OOOIDO 


CAU.  DQC32(VL.VH,DPlNDU,AilEA3) 

30  CPHAX  - 0.3DO'<AREAl«AR£A2-AReA3) 


SUBIIOUTINE  DCEtUXK.XXf  ,XXK,E) 
C PUSPOSE 


C CUL  DCELKXK.XXP.UK.E)  Co  osUuUcc  K and  E 
C CAU.  SCEU(XXP,XK,ICXXP,EP)  Co  csUulaca  K’  and  E- 
C DESCRIPTION  OF  PARAMETERS 


C METHOD 


INPUT  VALUE  OF  R'  IN  DOUBLE  PRECISION 
RESULT  VALUE  OF  K(k)  IN  DOUBLE  PRECISION 
RESULT  VALUE  OF  E(fc)  IN  DOUBLE  PRECISION 
RESULT  VALUE  OF  K'<k)  IN  DOUBLE  PRECISION 
RESULT  VALUE  OP  E'<k)  IN  DOUBU  PRECISION 


M.  Abraoowlcz  and  I.  A.  Stagun,  KANDBOOX  OF  MATHEMATICAL 
FUNCTIONS,  aevanch  prlnclng,  Uashlngcon  D.C.  : National  Bucaau 


PARAMnER<DPI-3.lalS92ES3SE979323SD0) 


COMFUTE  COMPLETE  ELLIPTIC  INTECRALS 
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C-(A-B)*.5D0 

BB-D5gilT(A*B) 


IF(DAB9(AA-B).LT.1.D-9*AA>  GOTO  2 


RETUKN 


SUBROUTINE  DCCWST(R,C1,C2) 


C CALL  DC0N9T(R,C1,C2) 

C DESCiaPTION  OF  FARAHETERS 


B - INFOT  VALUE  OF  R IN  DOUBLE  FRECISKW 

01  - BESULT  VALUE  OF  Cl  IN  DOUBLE  FRECISION 

02  - BESULT  VALUE  OF  C2  IN  DOUBLE  FRECISION 


C ROIABXS 


FUNOIION  FROGRAHS  REQUIRED 


C PREVIOUS  SUBROUTINE  AND  F\ 


uAlng  eqs.  (A-AB)  ami  (A-A9). 
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Q - -O.SDO*(BC  - DSqRT(BC*BC  «.DO*ACm;C]) 


SUNIOUTINB  DCF(AH,BH,CP) 


C CALL  SCF(AH,Bt(,CP) 

C DESCRIPTIQH  OF  FARAHETESS 


C AH  • ISPOT  VALUE  OF  «/h  IS  DOUBLE  PRECISIDR 
C BH  • INPUT  VALUE  OF  b/h  IB  DOUBLE  PSECtSION 
C CP  • RESULT  VALUE  OF  CP  IN  DOUBLE  PRECISION 
C REMARKS 


C SUBROUTINE  AND  FUNCTION  PROCRAHS  REQUIRED 
C PREVIOUS  SUBROUTINE  AND  niNCTIDN  PROCRAH  CALLS  REQUIRED 


IHPLICIT  REAL*B(A-H,D-E) 

COMXON  /MIST/  XX.XKP.XKK.XKKP.EK.EKP 
PARAMETER  (DFI-3.14159263358SVS32DO) 


THETA-ATAH(DEXP(W-At)) 

CSI^(BL00(l6.D0)-t4.D0»(AT+DU)GCDSlN<TOBTA))))/DPI 

•niA-DTAHH(AT) 


THE-&TAHH(»T) 

XlXniA*THa 

AA-l.M/eCOSH(AT) 

IF(DEX?(-1.01tlCWtI).NI.0.DCI)  Sa*l.DO/SCaSH<BT) 
XlP-BSQ!tT(XLP) 

CSl»2.M*DPI/(DLOe(l6.I>0)-2.»O*(BUM[THA)*DLOOaHB))) 

CSl'CSL/CHS 

CALL  DCEUCXL,XLP,XLL,E) 

CALL  [>CELI(i!LF.XL,XLL?,EP) 

CSl-2 -DO»XLL/XLLP 

CP-CSVCMS 

KmmN 


SUBROUTINE  DBU(E) 
C FUBPOSE 


DESCRimON  or  PARAHrrESS 

2 - INPUT  VALUE  OF  WPEOANCE  IN  DOUBLE  PRECISION 

XK  - RESULT  VALUE  OF  k IN  DOUBLE  PRECISION 

XKP  - RESULT  VALUE  OF  k'  IN  DOUBLE  PRECISION 

XKK  - RESULT  VALUE  OF  K(k>  IN  DOUBLE  PRECISION 

XKKP  - RESULT  VALUE  OF  K(V')  IN  DOUBLE  PRECISION 

EK  - RESULT  VALUE  OF  B(k>  IN  DOUBLE  FRBClSIffil 

EKF  - RESULT  VALUE  OF  E(k')  IN  DOUBLE  PRECISION 

DLNR  - RESULT  VALUE  OF  DLOC(k)  IN  DOUBU  PRECISION 
DLNXP  - RESULT  VALUE  OF  DLOC(k’)  IN  DOUBU  PRECISION 
REMARKS 

Villd  fsc  K/K'>0.  Oucpuc  ii  tranif«cred  vi«  comnn  block  /HIST/. 
SUBROUTINE  AND  FUNCTION  PROCRANS  REQUIRED 

PREVIOUS  SUBROUTINE  AND  FUNCTION  FROCRAH  CALLS  REQUIRED 


used  for  both  the  aodull  and  the  eoaplete  elllptio  Intestale. 
REFERENCES 

HATSEHATIK,  Vol.  7,  pp.  78-90,  1965. 


I,  HANDBOOK  OF  MATHEMATICAL 


P.  Byrd  and  H.  FrledBan,  HANDBOOK  OF  ELLIPTIC  INTECRAU 
FOR  ENGINEERS  AND  SCIENTISTS,  second  edition  1971.  pg.  298 


IMPLICIT  REAL'aCA-H.C-E) 


PARAHETER(DP1-3.1A1592653SB9793238DO,D1N2-.6931471805999‘>53DO) 

PA8AKETER(DP:OV2-0.5DQ*DPI) 
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PIC  - (Z.CT.rtA) 
C*RX 

IP(FIC)  0 - l.BO/RX 


give.  ch.  «.dulM  vUhln  0(l.D-l5)).  *o  th.  .arles  expsn.Ion! 


A-Ae2.DOeDASS(C] 

IP  (DASS(D).LT.DAJS(»)n.0-16)  GOTO  11 
1 COSTl(ra£ 

11  XK-(g/A)*(B/A) 

XXf-D5QRT(  1 . DO-XXeXK) 

C etlchaecu/geonetrlc  seen.  * 


DUIK-DLOC(XK) 

DU<KF>DLOG<XKP) 
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IF(DA8S<AA-B).LT.1.D-9»AA)  GOTO  22 

2 CONTINUE 
22  XKK-SPIOV2/A 

EK-(1.M  -.5DO*CS>*XKK 
XKKP-XKK/G 

EKP-DPI0V2/5KK-fJ[KKP*a.I10-BK/X») 


DU«-2.[10»61H2-DPI0V2/0 

DUKP-O.M 

nC-DEXP(DU«) 
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IP(FLS)  THEM 


DUKP-DIHK 


KEU.«g  FUNCTION  DFINBU(V) 
C PURPOSE 


C U-DPIHDU(V) 

C DESCRIPTION  OP  PAEANnERS 

C V . INPUT  VAUIE  OF  V IN  DOUBLE  PR3C1SKW 
C DPINDU  - RESULT  VALUE  OP  U(V)  IN  DOUBLE  PRECISION 
C REMARKS 


C SUBROUTINE  AND  FUNCTION  PROCRAMS  REQUIRED 

C PREVIOUS  SUBROUTINE  AND  FUNCTION  PROGRAM  CALLS 
C DELO,DCONST,DPARAN,DAREA.DQG32 
C METHOD 


REQUIRED 
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UPUCIT  HeAt'afA-H.Q-Z) 

INTECES  lER 

COHHON  /MIST/  XX,XKF,XKK,XKKP,EK,EKF 
COMMON  /AREA/  SU,CU,DU.SV,CV,DT,UL.UH,VL 


EU-BEL2<X  ,XKP , 1 . DO,  XK?*XKr) 
DFIHDV-02ERa(VL,VH,DH>FV,TOL,IER) 


ROUTINE  DFINDV’ 
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IKPUCIT  REALMS  (A-H.O-Z) 

COMHON/HSUB/KEPS , TARGET 
EXTERNAL  DVALV 


CALL  DELO(ZO) 


C l-DSqRT<  EKF/XXXP ) 

CALL  DPARAM(R,C1,C2.AP.BP,UD,A1I,BH) 


CALL  DqC32(0. DO, 1. DO, DVALV, AREA) 
EPre-CP*(RSPS-1.00)*AP*(l.D0-AP)/ 

& <l.DO-(l.DO-l.DO/REPS)MASEA-CP)/(l.DO-AP)) 
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HI-bA-  RATIO  OF  VIDE  CWDUCTOR  WIDTH  TO  SEPARATKHj 

VP  - PHASE  VELOCITY  REUTIVE  TO  THE  FREE-SPACE  VEIOCITY  OF 


IMPLICIT  REAL*B  (A-H,0-2) 

COMMON  /HIST/  XK,XXP,XKK,XKKP.EK,EKP 
PARAHETER(CC-2.«97924SED0) 

EXTERNAL  DVALV 


CALL  DB10{Z0)  ^ 


Cl-DSQRKEKP/XKKP) 


CAU  DPARAH(R,C1,C2.AP,BP.U0,AH,BH) 


CALL  [<1032(0.110,1. DO, DVALV, AREA) 
Cl-CP«(REPS-l.I>a)MP 

E (I.D0-(1-D0-1.D0/REPS)*(AREA-CP)/(1.D0 
EFFK-CI+CE 


)-AP)> 
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VP-1  .BO/DSCjRTlEFF*) 


SUBROUTINE  DPARAM(R,C1 


tip  parRM 


C SUBROUTINE  AND  FUNCTION  PROCRAHS  REqVIRED 
C PREVIOUS  SUBROUTINE  AND  FUNCTION  PSOORAK  CALLS  REQUIRED 


0 (A-IM)  with  B’-  b’/K'  ,A«  - «/h  using  (A-98),  BH  - b/b  using 
C (A- 101},  and  UO-uO/K  using  (A-U0>. 


C REFERENCE 


IMPLICIT  REAL»6(A-H,0-Z) 

COmON  /MIST/  XK.XKP.XKK.XKKP.EK.EKP 
PARAHETER[FIIHV  - 0.31S3098Siia379067D0) 


X - DSQBIT((C1*C1  • SXX)/(1.D0  - Cl*01))AX 
E - DEL2(X,XK,1.D0,SXX) 

F - DEL2(X,XK,1.D0,1.D0) 

SG  - DSORKCCl^Cl  • SXX)*(l.DO  - Cl»Cl)) 

SO  - SO«<1.DO/C1  • R/(l.DO  4-  C1*R)) 

AH  - (SO  + XETA)*J.D0*XKKP»P1I»V 


ELSEIFCR.BJ.l.DO)  THEN 
Call  R - l.DO,  B/K— >IHF1NITY 


X - DSQRI((e2*C2  - SXK)/(1.D0  - C2»C2))/XK 
E - DEL2{X,XX,1.1».SXK) 

5H  - DSQRT((C2*C2  • SXK)*(l.DO  - C2*C2)) 

BH  - (SH  - ZETA)*2.tlO»XIO!P*PIHJV 


R/DSQRT(1.I 
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SUKtOirriNE  BQC32  ie  froa  th«  IM  Systaii/360  Sctantlflc  Subrouclna 
Edlclen  (August  1970),  ps.30}.  No  DodifUatlon  of  thU  prograa  la  raqulted. 


SUBROUTINE  SSCD  (SN,CH,DH,X,CX) 
C PURPOSE 


C CALL  D5CS(SN,CN,l»t.X,aC> 

C DESCaiPIlOH  OP  PARAHETESS 

C SN  - RESULT  VALUE  OP  SK(X>  IN  DOUBU  PRECISION 

C CN  • RESULT  VALUE  OP  CH(X>  IN  DOUBLE  PRECISION 

C IH<  • RESULT  VALUE  OP  0K(X>  IN  DOUBLE  PRECISION 

C X - DOUBLE  PRECISKXt  ARCUHEHT  OF  JACOBIAN  ELLIPTIC  FUNCTIONS 

C CX  - CMPLEKENTART  MODULUS  IN  DOUBLE  PRECISION 
C REMARKS 


C PREVIOUS  SUBROUTINE 


FUNCTION  FROSBAH  CALLS  REQUIRED 


C REFQIENCE 


0 R.  Sullrseh,  NUMERICAL  CALCULATION  OP  THE  ELLIPTIC  INTECRAIS 
0 AMD  ELLIPTIC  FUNCTIONS. 

C HANDBOOK  SERIES  OF  SPECIAL  PUNCTKWS 
C NUMERISCHE  HATHEHATIK  VOL.  7,  19(3,  PP,  78-90, 


DIMENSION  EM(13), 01(13} 


INTECEH  l.II,L 


iF<B<c.tie:.o.iio>‘niEN 

FliS-(EHC.LT.O-SO) 


IVDSQKTCfi) 


EKODSQaKEMC) 

cc-o.5i)o*(Ma«:) 

1F(DABS(A-EMC).LE.1.0- 


SN-SSIH<U) 

CN-DC03<U) 

A-CN/SN 


I>H>(EN(II]tA)/(BM) 

CONTINUE 

A-l.B0/IsqBT(C»*2*l  .M) 


1F(SN.LT.0.D0)THEN 


CN-1.D0/DCQSH(U) 

SH>DTiUn«U) 


{k  ne.r  l.M). 


HEAl^S  FUKCniOR  DVALV(D) 
C PUfIFOSe 


C V-0VALV(U) 

C SESCUniON  Of  fAKAHETESS 

C U - DOUBLE  PRECISION  VALUE  OF  U-u/K' 

C V - RBSULTWC  VALUE  OF  V-v/K'  IN  DOUBLE  PRECISION 

C BSHARKS 

C SUBBOUTIME  AND  PUHCTION  FROGRAHS  REgUIRED 
C PREVIOUS  SUBROUTINE  AND  FUNCTION  PROCRAM  CALLS  REQUIRED 
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IKPUCII  MAL*8  (A-H.O-Z) 

COHKCM  /HIST/  >X,XKP,XKK,XKK?,EK,EK? 

CALL  DSCD(SU.CU,I>U,V»XKK,XX?) 

E-DEL2  <X , XKP , 1 . M , XKP* W ) 

FACT-.5K>*3.U15926535a97932DO/W<KP»ll.DO-U)tU*EK-E 

X-OSQRI(FACT/{SU»{CU»Di;-Sl)»PACT)))/Tm 

F-DEL2(X.XK,l.DO,l.DO> 


REAL'S  FUNCTION  DXOFU(U) 
C PURPOSE 


C P(U)  - nXOFU(U) 

C DESCRIPTION  OP  PARAHRTCKS 

C U • INPUT  VALOR  OF  U IS  DOUBLE  PRECISION 

C DXOFU  - RESULT  VALUE  OF  F(U)  IN  DOUBLE  PRECISION 


C SUBROUTINE  AND  FUNCTION  PROCRAMS  REQUIRED 


C PREVIOUS  SUBROUTINE  AND  FUNCTION  PSOCRAK  CALLS  REQUIRED 


IMPLICIT  REAI;*8(A-H,< 
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COMMON  /MIST/  XK,XKP.XKK,XKK?,EK.,EKP 
COMMON  /»R£A/  SU,CU,DU,SV,CV,DV,UL,IIH,VL, 


BDF-  T3»SV2*CV2»(CU*Cn  -l-  D1)*DU*DV2)/T2 
AAK-  *»CU*BU»m*CCV2  - Tl) 

BEN-  (R*SU»DVl  ♦ I2)»»2  t (»»CU«BU»5V*Cn)»»2 
XH-  (S**2»<ACP-8DF)4AAK)/DB« 

EF-  I3»S»2/T2 

e-DEU  ( X , XXP , 1 . DO , XKP*XKF) 

XK-  U ♦ 2.0O»XXKP».31g3O98a618379067DO»(E-El!*tttXHtEF) 


R£AL*8  FUNCTION  DXOFV(V) 
C FUBPOSE 


F(V)  - DXOFV(V) 
OESCRiniON  OF  PARaHETEBS 


N DOUtU  PRECISION 


C SUSSOUTINE  AND  FUNCTION  PBOCRAMS  REQUIRES 
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C PREVIOUS  SUBROUTIKE  AKD  nmOTION  PROCRAH  ' 


REgUlRED 


C REFERBICE 


IKPUCIT  REAL*8(A-H.Q-Z) 


CAU.  nSCO(SV,CV,DV,V*X)CKP,XK) 

CTl-CV/XK 

BVl-DV/MC 


II-  SU*SU»SV2 


IF(R.NZ.O.DO)  THEN 


ACF-  T3»DV2*(CV2  - Tl)/T2 

BDF-  T3»SV2»CV2»(CU»CU  * DU*DU»DV2)/T2 

AAK-  R*CUi»IIU»DVl»<CV2  • II) 

DEN-  (R»SU*DV1  ♦ T2)«2  * (»»CU»DU*SV*CV1)»*2 
XH-  (R**2»(ACF-BDP)tAAX)/DEN 

BF-  T3*SV2/I2 

XK-  U t 2.»O*XKKP».3U3O988Sia379067DI)»(EU-EK*U*XH»EP) 
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RZU.*I  FUNCTION  DZERO<AX,BX,F,TOL,ZER} 


IF(I>SICM(1.IH),FA)'>'SSIGN(1.I 


DO.FB))  100,20,100 
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